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FOREWORD 


This report is the second of two irtiich have resulted from contract 
NAS8-33097 between Physical Sciences Inc. and the NASA Marshall Spaceflight 
Center, for a research study entitled "Laser-Heated Thruster." The first 
report is an Interim Report, dated February 1980, PSI TR-205. The present 
report is the Final Report PSI TR-220. 

The effort accomplished under the contract, over the period 
September 1, 1978 through September 30, 1980, is documented in these two 
reports. Together they present the technical work performed. The Interim 
Report contains most of the work on laser absorption physics, radiation 
physics, thermodynamics and transport properties, and the preliminary ver- 
sion of the thrust chamber design computer program. In this Final Report 
the first five chapters contain a summary of all the work, which can serve 
as an Executive Summary. The nine appendices of this Final Report contain 
the details of additional technical work not presented in the Interim Report. 
The two reports should bo read together to obtain complete documentation 
of the work accomplished. 
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I. INTRODUCTION 


This report is the fifth in a series of reports vxitten for NASA 
by Physical Sciences Inc. on the subject of CW laser propulsion cuid laser- 
heated thrusters. The first three reports were written under sponsorship 
of the NASA Lewis Research Center, while the fourth and the present (fifth) 
ones were written under sponsorship of NASA Marshall Spaceflight Center. 

The first report (Ref. 1) presented a study of the mechanisms by 
which laser energy can be absorbed into a gas. It also reported on a study 
of the flow of an ideal gas in a nozzle, when undergoing laser energy absorp- 
tion. The stability of this absorption process to small disturbances was 
also considered. 

The second report (Ref. 2) dealt primarily with the flow of hydrogen 
which IS adjsorbing 10.6 pm laser radiation by inverse Bremsstrahlung. The 
hydrogen was treated as a real gas in chemical equilibrium. The flow was 
quasi-one-dimensional. Radiation losses from the hot hydrogen were included 
by studying its radiative properties and constructing a sinple model to use 
in the flow equations. The pure hydrogen laser-heated thruster requires 
the presence of a laser supported combustion (LSC) wave at the entrance to 
the heating chamber to provide the absorption mechanism. In Ref. 2 the prop- 
erties of this wave were estimated rather crudely. 

A third report (Ref. 3) focused on the structure of LSC waves in 
equilibrium hydrogen absorbing 10.6 ym radiation, to provide more detailed 
information on this important element of pure hydrogen thrusters. A one- 
dimensional model was constructed, including radiation losses and axial radia- 
tion transport. Several two-dimensional effects were incorporated into this 
one-dimensional model, including laser beam convergence and transverse con- 
duction and radiation losses. The relation between laser intensity and mass 
flow per unit area into the wave was determined numerically. A one -dimensional 
expansion of the heated hydrogen to a nozzle throat was also calculated for 
one case. Peak hydrogen temperatures of 20,000 K were found in LSC waves. 
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These studies for NASA/Lewis led to the general conclusion that laser 
energy could be used to heat hydrogen in a stable mamner, but the resulting 
gas was very hot and would lose a substantial portion of its energy by radia- 
tion. It became clear that further studies should focus on ways to reduce 
the heat losses and the associated wall-cooling problems. 

At the same tune, the sponsorship of laser propulsion work within 
NASA was shifted to the Marshall Spaceflight Center, under whose aegis PSI 
has continued their studies. The present report is the second under Marshall 
sponsorship, representing the final report on a two-year research study. 

This work has explored some alternatives to the pure hydrogen rocket, which 
will permit operation at lower temperature with smaller losses and less- 
stringent wall-cooling requirements. It has developed two-dimensional (axi- 
symmetric) computer programs useful for the design of the thrust chamber 
of a laser-heated thruster. 

One approach to reducing the gas temperature is to seed the hydrogen 
with a small amount of easily lonizable gas, such as an alkali metal-vapor. 
This seed will produce electrons at a much lower ten?>erature than hydrogen- 
Since electrons are a necessary element for laser energy absorption by inverse 
Bremsstrahlung, this will permit ^lbsorptlon to start at a lower temperature 
than in pure hydrogen, and so result in a lower tenperature gas. An LSC 
wave IS still needed to heat the gas from the inlet tenperature up to the 
teirperacure at which sufficient ionization has occurred to permit laser energy 
absorption. This LSC wave operates in seeded hydrogen, and will have a dif- 
ferent wave speed vs. intensity relationship than a pure hydrogen LSC wave. 

One of the tasks in the present study was modification of the pure hydrogen 
LSC wave model developed in Ref. 3 to include the presence of cesium vapor 
seed. 

A second approach to reducing the gas tenperature is to seek laser 
energy absorbers which operate at the inlet gas tenperature and continue 
absorbing all the laser energy as the gas tenperature rises. This requires 
a low tenperature absorber as well as a medium temperature absorber. The 
latter cam again be an alkali metal vapor, such as cesium. The former must 
be a molecular absorber "tuned" to the particular laser radiation being used. 



Candidate molecular absorbers were identified based on the studies of Ref. 1. 

A likely candidate for low temperature absorption of 10. C im radiation is 
water vapor, which will absorb this radiation until it dissociates around 
4,000 K, after which the alkali metal vapor will produce enough electrons 
to let inverse Bremsstrahlung absorption become effective. With both low 
and medium temperature absorbers, there is no need for the use of the LSC 
wave mechanism to provide the initial gas heating. The absorption proceeds 
continuously from the inlet temperature until all the laser energy is absorbed. 
This heating mode will be referred to as continuous heating, in contrast to 
the LSC wave heating mode, where only medium and high tenperatura (hydrogen) 
absorbers are used. 

The continuous heating mode can also be used with other laser wave- 
lengths. For CO lasers, with radiation around 5.3 )jm. Ref. 1 suggests that 
either NO or H^O can function as the low temperature absorbers. In this 
case,- there is also a molecular cibsorber, CO, which survives to medium tenpera- 
ture (at least 6,000 K) , and could be used instead of the al):ali metal vapor 
seed. It is advantageous to avoid this latter seed, if possible, because 
Its high molecular weight reduces the rocket specific inpulse. 

Introducing seeds into the primary propellant gas conplicates the 
gas storage and supply system, since it requires going from a mono-propellant 
to a bi-prcpellant or tri-propellant system, although the seeds do not have 
to be kept separate from the primary propellant in this case. There is another 
possible cipproach to reducing the temperature which does not require any seeds. 
This idea is to provide a co-flowing buffer gas around the central laser- 
heated hydrogen core, which can act as .an insulator for the walls. After 
the core is heated, it will mix with the buffer and finally provide a uni- 
formly heated gas flow at a lower temperature to be expanded and provide 
thrust. The buffer gas can also be hydrogen at the same inlet conditions 
as the core gas. The only difference is that the laser beam diameter is 
smaller th.an the gas stream diameter, so only the central part of the gas 
IS heated by the laser. If this approach yields satisfactory specific im- 
pulse, and low enough wall heating to be handled, it is perhaps the simplest 
emd most desirable one . 
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The present study was undertaken to provide the tools for assessing 
specific examples of these three alternate approaches to CW laser-heated 
thrusters. We h.ive considered the following combinations of laser wavelength 
and propellant gas mixtures; 

1. CO^ (10.6 '.ai) radiation with H^/Cs (LSC wave heating mode) 

2. CO^ (10.6 tm) radiation with H^/H^O/Cs (continuous heating 

mode) 

3. CO (5.3 Un) radiation with H^/CO/NO or H^/CO/H^O (continuous 
heating mode) 

4. CO^ (10.6 .an) radiation with (LSC wave heating 

mode) 

To evaluate each of these altcrnati\-es the following physical prop- 
erties of the gas, mixture must be modeled: 

1. Laser energy absorption properties. 

2. Radiative and absorptive properties. 

3. Thermodynamic and transport properties. 

These gas properties must be incorporated in a gas flow model whj.ch 
includes the relevant effects of heating by laser energy absorption, gas 
radiation emission and absorption, and the usual flow effects of ccnvec- 
ticn and heat transfer. 

In the interim report issued under the present study (Ref. 4) most 
of the physical properties of the gas mixtures were presented. The absorp- 
tion of 10.6 and 5.3 ..m radiation by inverse Bremsstrahlung, and by water 
vapor, and the absorption of 5.3 tm radiation oy H.,0, CO and N'O were described. 
The radiati%’e and absorptive properties of the important radiating spocies in 
the gas mixtures were f^cdelod. The thernody.namic and transport properties 
of a mixture of hydrogen and cesium up to 20,000 K '^ere developed, for use 
in LSC wave calculations. (For the other gases, in the continuous heating 
mode, 3ta.ndard chemical equilibrium codes can be used, sinco the temperatures 
of interest for ccnti.nuous ■'eating do not go -uch beyond 5, COO K. ) 



Also contained in the interim report was an ax i symmetric flow model 
for the continuous heating mode, including a laser beam shape, and a radiation 
model for the particular case of a 30 atm mixture of H 2 /H 2 O/CS = 0.945/0.05/ 
0.005. 

The present final report presents additional work performed since 
August 1979. It should be looked upon as supplementary to the interim report 
(Ref. 4) , and does not repeat the material of that report. The two reports, 
taken together, present the work accomplished under the present contract 
with NASA/MSFC. 

Chapter II of this report describes some results obtained using simple 
flow models for the LSC wave heating cind continuous heating modes. These 
models can be used for mapping out the regions of interesting thrusters and 
providing order-of-magnitude estimates of their performance and heating prop- 
erties. They provide very useful simple tools for parametric e valuation of 
laser-heated thrusters. The details of the models are given in Appendices 
A, B, C and D. 

Chapter III describes the axisymmetric flow thrust chamber design 
codes for both heating modes, as well as the LSC wave code whose output is 
needed for the LSC wave heating mode design code. The details of the flow 
models are provided ir Appendices E, F, G, and H. 

During the early parr of the present study, a literature survey was 
made of the eibsorption coefficient of water vapor for 10.6 Um radiation. 

It was found that the literature showed a steep minimim near 500 K, predicted 
theoretically, and a number of emission experiments wh’ch tended to confirm 
this minimum. However, there appeared to be no laser absorption experiments 
on water cibove room temperature. Then a report appeared (Ref. 5) in which 
laser absorption measurements were made in a laser-heated gas. The values 
of absorption coefficient presented there were an order of magnitude or more 
above the values in the rest of the literature. In view of the inportance 
of water vapor as a possible low ten^jerature absorber, it was decided by 
NASA/MSFC and PSI to perform a laser cibsorption experiment on water vapor 
heated in a shock tube , where the gas conditions could be predicted very 


reliably. The results of this experiment are presented in Chapter IV of 
this report; the details of the experiment are given in Appendix I. 

The conclusions which can be drawn from the work performed during 
the present study are set forth in Chapter V» 


II. SIMPLIFIED MODELS OF LASER-HEATED THRUSTERS 


Parametric studies of laser-heated thrusters can be facilitied by 
the use of simplified models of the thrust chcunber, which contain the essen- 
tial physical elements of the thruster, but can be in^lemented in small, fast 
computer programs. Such models can be used to map out the operating charac- 
teristics of laser-heated thrusters, and describe approximately the regions 
of power, performance and wall heating in which they operate. These maps 
can then be used to indicate how engineeri.ig constraints, such as maximum 
wall heating rates, limit the performance of thrusters which can be built. 
This, in turn, provides information about what thrusters are likely to lead 
to useful rocket engines. The maps can also be helpful in designing experi- 
ments to prove the feasibility of laser-heated thrusters. 

In order to provide this important information, PSI has developed 
simplified models for both the LSC wave and continuous modes of gas heating. 
These models are described in Appendices B and A, respectively. They relate 
the following thruster parameters; 

F = Thrust 

I = Vacuum specific impulse 
sp 

P = Incident laser power 

Lj 

= Average wall heating rate 

R = Thrust chamber radius 

L = Thrust chamber length 

The models are implemented in simple, fast computer programs which 
use the thermodynamic state of the gas and some of the above parameters as 
inputs. They use conservation of mass, momentum and energy to calculate the 
other parameters. Included are estimates for convective and radiative wall 
heating. The thrust chamber geometry is modeled by a cylinder. In the con- 
tinuous heating mode the wall heating is calculated by taking the gas tempera- 
ture to be uniform at the highest temperature, which is the exit temperature. 
In the LSC wave heating mode, the losses from the LSC wave are also calculated 



using the highest temperature in the wave, and those from the mixed gas 
region are found using the exit ten^serature . These are believed to be 
conservative estimates of the losses. On the other hand, the wall heating 
rate is the average over the chamber walls, and does not consider the pos- 
sibility of hot spots along the chamber length. 

The results of these models are expected to be accurate to better 
than am order of magnitude, and perhaps as good as a factor of two or three. 
With that understanding, they should be very useful in mapping out thruster 
parameters. 

The models have been exercised to provide some examples of the re- 
sults which can be obtained. Figure 1 shows a map for an buffer 

thruster absorbing 10.6 ym radiation in the LSC wave heating mode at a 
pressure of 3 atm. It covers the range of 5 kN to 0.5 kN thrust (1124 to 
112 pounds), and of 1094, 1402 and 2007. The combination of laser in- 
tensity I and mass flow per unit area m is appropriate to a hydrogen LSC 
X> A 

wave, as calculated in Ref. 3. The laser power required varies from 3 to 

2 

60 MW, cmd the average wall heating rate from 0.8 to 4 HW/m (80 to 400 
2 

W/cm ) . If one recalls that the power in the exhaust is P„ = F I /2, 

G g sp ' 

then It varies from 2.7 MW at 0.5 kN, 1100 s to 49 MW at 5 kN, 2000 s. 

Thus a large percent of the laser energy goes into the gas, making for an 
efficient absorption process. The heating rates appear to be in a range 
that can be handled with fairly conventional cooling methods. The wall 
heating is dominated by radiative rather than convective contributions. 

A similar nap for 30 atm is shown in Fig. 2, with the appropriate 
mass flow for the same laser intensity. It shows very similar results, 
with some reduction of heating rates at the same thrust and primarily 

because the LSC wave is thinner at the higher pressure, because of enhanced 
laser energy absorption coefficients. The decrease in LSC wave size also 
changes the relative contributions of radiation and convection to the wall 
heating, making the convective contribution dominate. 

Turning to the continuous heating mode, we have calculated some 
cases for a gas mixture of H^/H^O/Cs = 0.945/0.05/0.005 absorbing 10.6 ym 
radiation. This mode does not have a relation between intensity and mass 


8 









3157 



Qw(W/m^) 


Figure 2 Laser power and average heating rate nap for 
H^/H^ I^C wave heated thruster at 30 atm. 
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flow per unit area, so instead the radius of the thrust chamber is chosen. 

Figure 3 is for 3 atm, R = 0.1 m, for F = 0.1 to 5 kN and = 1042, 1390 

and 2040 s. This map falls in a similar power and heating range to Figs. 1 

2 

amd 2, varying from 0.8 to 50 MW power and 0.7 to 8 MW/m heating rate. At 
the same thrust and specific impulse, this case has in fact somewhat higher 
heating rates at the higher thrusts than shown in Fig. 1. The heating here 
is mostly dominated by convection, except at very low values of thrust. 

Raising the pressure level to 30 atm gives the map in Fig. 4, which 
is shifted almost a decade to the right compared to the previous maps, being 
centered on 10 MW/m instead of 1 MW/m . However, these higher heating rates 
should be viewed with some caution. At 30 atm the length of these thrusters 
is very small, varying from 2.6 cm at 3500 K to 0.6 cm at 6000 K. These small 
lengths are caused by the large absorption coefficients at 30 atm, since the 
length is taken as the inverse of the absorption coefficient. They in turn, 
lead to average high convective heating, which is inversely proportional to 
the square root of the length. Longer thrust chambers might have smaller 
average heating, although the radiative component of the heating will also 
change. The laser power levels of these thrusters are very much the same as 
those at 3 atm in Fig. 3. The heating is a small fraction of the laser power, 
so considerable changes in heating do not alter the power level much. 

These calculations are meant to illustrate the kind of results which 
can be obtained rapidly with the simplified thrust chamber models. They are 
by no means a parcimetric study. They do show what values of laser j>ower are 
needed, and what average heating rates have to be handled for thrusters which 
might be of operational size. The powers are not much larger than the exhaust 
gas power, and the heating rates seem within a range which are dealt with in 
present rocket practice. 

Another use for the simplified models is in the design of small-scale 
experiments. As an example of this sort of calculation, we have considered 
thrusters which can be designed with a 30 kw laser beam. Again, calculations 
have been made for both heating modes at 3 and 30 atm pressure. 

Table 1 shows the results for the LSC wave heated mode at 3 atm, with 
9 2 

laser intensity of 3 x 10 W/m , ^lnd I = 1094, 1833 and 2454 s. The thniat 

sp 
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levels are a few newtons, cind the average wall heating rates are 2.45 to 
2 2 

13.8 MW/m (245 to 1300 W/cm ) . Also shown in Table 1 are the thruster 

radius R^, the laser spot size R^, the ratio of power going into the wall 

to that in the gas at the thrust chamber exit P /P_ and the ratio of radia- 

w G 

tive to convective wall heating Q /Q . It can be seen that these small size 

R C 

thrusters are not as efficient in putting power into the gas as are the large 
ones previously considered, and they receive about equal amounts of radiative 
and convective wall heating. But the size and heating rates are in the range 
of possible experimental hardware. 


A similar calculation for 30 atm is shown in Table 2, for the same 

30 kW, but with the lower intensity LSC wave that is permitted at 30 atm, 

9 2 

namely 10 W/m . The over-all radius is about the same, but the spot size is 
larger because the intensity is lower. The power efficiency J-s also 

better, and the wall heating is convection dominated, because the shorter 
lengths associated with the higher absorption coefficients at 30 atm result 
in a smaller radiating volume. 


The small-scale continuous heating node is represented by the re- 
sults shown in Figs. 5 and 6. They are for the same gas mixture as before, 
Ht/H 20 /Cs = 0.945/0.05/0.005, and for a thruster radius of 2 cm. Figure 5 
presents the 3 atm case down to F = IN. The 30 kW power level will pro- 
duce about 2N at 1042 s or IN at 1389 s. The corresponding average wall 
heating rates are 0.3 and 0.5 MW/n“, wnich are quite low, although the 
thrust chambers are between 37 and 55 cn long at these low temperatures. 

2 

At 30 atm. Fig. 6 shows that 30 kW will produce 5N at 882 s for 1.7 MW/m , 

2 

oi 2N at 1570 s for 9 >W/m . At this pressure, the thrust chambers are 2.6 
and 0.6 cm long, respectively. The results indicate the possibility of pei- 
forming small-scale, continuous heating thruster experiments within reason- 
able engineering constraints. 

This small group of examples shows the usefulness of the smylified 
thruater models as engineering tools for defining both full-size and small- 
scale thruster parameters. Furthermore, tne calculations show promise that 
laser heated thrusters with measuralple thrust can be tested on a small scale, 
and that large-scale thrusters can be made with useful values of thrust and 
specific impulses over 1000 s. 
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III. THRUST CHAMBER DESIGN PROGRAMS 


One of the objectives of the present study was to develop computer 
programs which could be used for the design of laser heated thrust chcimbers. 
Because of the presence of the LSC wave, the program for the LSC wave heat- 
ing mode differs from that for the continuous heating mode, so that two pro- 
grams are needed. In addition, a program for calculating LSC wave structure 
IS needed, since the output of that program is used as the input for the 
thrust chamber design program in the LSC wave heating mode. 

The thrust chamber design programs solve the axisymmetric flow of 
a gas in chemical equilibrium, which is undergoing laser heating, radiation 
energy loss and radial transport, and conductive radial transport. The 
equations of motion are written using the boundary layer approximation that 
the axial component of conductive energy and motrentum transport can be ig- 
nored compared to the radial component. The equations use the stream func- 
tion as the independent variable in place of the radius. Only the initial 
radius of the thrust chamber is given at x=0. The shape of the chamber, 

R (x) , IS determined as a result of the calculation, and is mainly influ- 
w 

enced by the axial pressure gradient, which is input as a function of axial 
distance. 

In addition to the chamber radius R^(x), the programs calculate the 
local heating rates at the walls, the velocity, temperature and laser in- 
tensity profiles, and the total energy transferred to the walls up to any 
axial station. 

The equations for the continuous heating program (called LHTE) are 
presented in Appendices E and F. The equations for the LSC wave heating 
program (called CLHTE) are the same as those for LHTE, but the initial con- 
ditions at x=0 are different. These initial conditions are obtained from 
a solution of the LSC wave which stands at the entrance of the thrust chamber. 
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The program which solves for the LSC wave structure in hydrogen, or 
a mixture of cesium and hydrogen, is called LSCWCS. The equations which 
describe it are presented in Appendix G. They are one-dimensional flow 
equations for a constant area flow of a laser heated gas in chemi''al equi- 
librium, with radiation loss and aicial transport, and conductive axial 
transport. They include laser beam convergence effects and esiimates for 
radial energy loss. 

The method of coupling the output of the LSC wave program LSCWCS 
to the thrust chamber code CLHTE is described in Appendix H. The gas em- 
anating from the LSC wave at some chosen axial station in the wave (currently 
the maximum temperature point) is used as the initial state for the central 
core of the thrust chcmber. An annular region around this core is taken to 
be a buffer region of gas which is not heated by the laser. The initial tem- 
perature of this buffer gas is higher than the inlet temperature of the gas 
entering the LSC wave, because the radial energy losses from the LSC wave are 
used to heat up the buffer gas as it flows past the LSC wave. This slightly 
heated buffer gas is the initial state for the outer cuinulus of the thrust 
chcimber, at the same axial station as the LSC wave heated gas. As the two 
gas streams flow through the thrust chamber they mix to some temperature in- 
termediate between the buffer gas inlet temperature and the high temperature 
of the laser heated gas. 

A few sample calculations have been made with these codes, and some 

results will be presented. First consider the LSC wave heating mode, using 

5 2 

hydrogen and a hydrogen buffer, for 3 atm, with a laser intensity of 10 W/cm . 

The LSCWCS code yields a gas at about 17,400 K, and a residual laser inten- 
4 ? 

sity of about 2 x 10 W/cm , so 80% of the laser energy is absorbed in the 
LSC wave up to the peak temperature station. The laser spot radius is 
3.99 cn, and the initial chamber radius is 8.2 cm, so the buffer gas is 
contained between 3.99 and 8.2 cm. The initial buffer gas temperature is 
calculated to be 433 K, representing a 133 K temperature rise from the inlet 
temperature of 300 K, caused by heating from the edges of the LSC wave. 
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The CLIITE program was run for 12.26 cm in the axial direction. Radial 

profiles of the laser intensity are presented in Fig. 7 at the three axial 

stations x = 0, 0.305 and 12.264 cm. The initial profile shows the initial 

intensity coming through the LSC wave up to a radius of 4 cm, and no laser 

in the buffer gas region. At 0.305 cm, the intensity has dropped to about 
2 

5700 W/cm in the core, as further absorption takes place. (The more gradual 
transition to zero intensity at this station is an artifact of the numerical 
solution.) By 12.264 cm, the laser intensity has vanished, all the energy 
having been absorbed. 

The corresponding temperature profiles .ire given in Fig. 8. The 
initial profile shows the discontinuous profile from 17,400 K in the core 
to 433 K in the buffer region. At 0.305 cm, a slight further heating of the 
core has occurred, from the residual laser intensity at the back of the LSC 
wave, and the core temperature has risen to 19,000 K, with very little dif- 
fusion of heat to the butfer gas as yet. By 12.264 cm cooling and diffusion 
have begun to show their effects, and the temperature profile has begun to 
smooth out, with its naximiun down to 15,500 K. 

The axial velocity profiles arc shown in Fig. 9. They show little 
change between 0 and 3 cm, both giving nearly discontinuous profiles from 
the 30,000 cn/s emanating from the l£C w.ive to the 235 cm/s at which the 
slightly heated buffci gas flows. (The inlet gas is flowing at U'2 cm/s.) 

By 12.264 cm the velocity profile has accelerated in the center and diffused 
outward, with a maximum velocity of 36,000 cm/s. 

A plot of the wall radius vs axial distance is shown in Fig. 10. 

Recall that this is calculated from the input prcssuie gradient distribu- 
tion, which was taken as constant tot this case at dp/dx = -2.3 \ 10 ^ atn/cn, 
resulting in a nc<]l igiblo-'change in pressure up to 12.264 cm. The wall pro- 
file shows a decrease from 8..’ cm at the inlet to 5.2 cm at 12 cn, an area 
docioaso from 211 to 85 cm", or a factor of 2.5. The profile looks quite 
smooth, except for a slight bump near 0.8 cn. 
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Figure 7 Radial profiles of laser intensity for LSC wave 
heated thruster from program CLHTE. 



Figure 8 Radial profiles of temperature for LSC wave 
heated thruster from program CUHE. 
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Fiquro 9 Radial profiles of axial velocity for LSC wave 
heated thruster from program CLHTE. 
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The power in the vjns as a function of axial distance is shown in 

Fig. 11. It starts at the value emanating from che LSC wave of 3.46 MW, 

rises slightly to 4.3 MW at x = 0.6 cm while the rest of the laser energy 

IS being absorbed, and then decays slowly as it loses energy to the walls. 

It reaches 2.59 MW at 12.264 cm. The laser power input to the LSC wave 

5 2 

in a spot of radius 3.90 cm with an intensity of 10 W/cm was 5 MW, so 

about 90% came through the LSC wave in the gas or in the residual laser 

beam, and about half is still in the gas at the last station. However, 

for this sample case, the gas is not yet accelerated to the throat, since 

the sound speed at 3 atm, 15,500 K for hydrogen is about 1.4 x 10^ cm/s, 

4 

while the maximum gas speed is only 3.6 x 10 cm/s. To make this case more 
efficient, it would be necessary to accelerate the gas much faster, which 
would require a larger input pressure gradient. 

Finally, the distribution of wall heating rate is shown in Fig. 12. 

3 2 3 

The peak value is about 6 \ 10 W/cm at 1 cm, with a slow decay to 2 x 10 

W-'cm“ at the last station. This cctresponds to about 1.99 MW put into the 

walls. The oneigy balance has then 4.57 MW in the gas or walls at the last 

station, while the jxjwer emanating from the LSC wave was 4.46 MW. The s.ight 

discrepancy of 2.5% repiesents the numerical inaccuracy in doing the various 

axial and radial infograls to find the integrated power contributions. 

The results of an example using the continuous heating proaram LUTE 
are shown in Figs. 13-19. This case is for a mixture of H.,,H^0/Cs = 
0.945/0.05 0.005 at '0 atm entering at 1000 K, with a laser jx^wer of 1 MIV 
incident on a thrust chambet whoso initial radius is 12.39 cm. The shape 
of the inccning beam is rel.ited to the Hessel function, as described in 
Eqs. (5-12) and (5-13) of Ref. 4, with a {-’oak intensity at the center line 
of 3100 W'en". This calculation has been run to 20 cn. The figures have 
profiles at 0, 8.971, 10.231 and 16.482 cm. 

Figure 13 gives the laser profiles. Tlie initial shape is shown at 
X = 0, and the subsequent profiles show the laser energy being absorbed well 
in the center of the chamber, but not as well near the walls, where the tem- 
perature IS cooler because the walls are held at a fixed temperature of 1000 
K. Finally, at 16.432 cm the laser energy is almost completely absorbed. 


even near the wall. 
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Figure 13 Radial profiles of laser intensity for continuously 
heated thruster from program LHTE. 
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Figure 14 Radial profiles of temperature for continuously 
heate^i thruster from program LHTE. 
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Fiqure 15 Radial profiles of axial velocity for continuously 
heated thrustci froin proaraia LHTE. 
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The temperature profiles are presented in Fig. 14, and show an in- 
crease from the 1000 K inlet to 4000 K, and then a reduction to 3400 K as 
the energy transfer to the wall continues to take energy from the gas after 
the lasct energy is fully absorbed. 

The velocity profiles of Fig. 15 show a steady increase from the 
inlet velocity of 9.3 cm/s to a value of 150 cm/s at 16.482 cm. Of course, 
this factor of 15 acceleration still represents only a small velocity, and 
the reason can be sought in the wall radius distribution given in Fig. 16. 

It decreases rapidly from the initial value of 12.39 cm to 7 cm at x = 2.5 cm, 
and then rises to a peak after a further shallow minimum. Following the 
peak It decreases again. This shape is a result of the pressure gradient 
distribution input, which is dp/dx = -2.745 x 10 atm/o-, a constant. It 
IS clearly not a realistic thrust chamber shape. The program must be 
exercised with different input functions dp/dx to learn how to make reason- 
able nozzle shapes. 

The distribution of power in the g. s is shown in Fig. 17. It rises 
as the laser energy is absorbed, and then decays as losses to the wall con- 
tinue. It IS 0.172 MW at 16.482 cm. The wall heat rate distribution is 
given in Fig. 18. It rises rapidly at the beginning because the rapid de- 
crease in wall radius causes unabsorbed laser energy to hit the wall, which 
IS by far the largest contribution to the wall heating there. The slow de- 
crease after x = 0.7 cm is caused by the change i.i wall angle shown in 
Fig. 16, which means a smallei fraction of the laser intensity is normal 
to Che wall. The rapid decrease starting at 3 cm again reflects a smaller 
wall slope seen in Fig. 16. When the wall radius begins to increase again, 
there is no longer laser heatirg of the wall since the beam is now blocked. 
However, the wall heatJng from the gas is beginning to rise as the gas heats 
up. This accounts for the increase in heating up to about 11 cm, where the 
wall radius peaks. From there on there is a slow further rise, and then 
a decay as the gas cools. 
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At 16.482 cm, there is 0.172 MW in the gas, 0.004 MW left in the 
laser beam, and 0.984 MW has been put into the wall. The total is 1.17 MW, 
compared to the input energy of 0.993 MW in the laser beam and 0.27 MW in 
the incoming gas flow, for a total of 1.02 MW incoming. This 15% discrep- 
ancy in energy balance is again caused mostly by inaccuracies in evaluating 
the integrals. However, it should be noted that the program does not con- 
template walls which increase in radius, and the calculation of the laser 
energy being put into the wall is no longer correct when the wall radius 
begins to increase. 

These two examples show the kind of results which can be obtained 
with the LHTE and ClilTE programs. However, the programs have not yet been 
exercised extensively enough to provide the experience necessary to make 
them useful tools for thrust chamber design. The relation between the input 
pressure gradient distribution and the output wall radius distribution must 
be explored, for example. When a body of computational results has been 
obtained and analyzed, these programs should provide the information neces- 
sary to continue the development of laser-heated thrusters. 

A useful way to run those programs would be in conjunction with the 
simplified models discussed in the previous chapter. Those models could be 
used to -define particular thrusters which seem interesting, and then the 
more elaborate programs described in this chapter could be used to provide 
the detailed information on those thrusters. 

It should be pointed out that some of the modeling included in the 
three programs LSCWCS, LHTE and CLHTE needs improvement. One such area is 
the radiation modeling for H^/Cs mixtures in LSOTCS and CLHTE. The present 
radiation model in the LSC wave program does not include cesium radiation, 
which may be an important radiator. In the thrust chamber program which 
couples to LSCl^CS, namely CLHTE, the radiation model only goes up to about 
6000 K. The effect of high temperature radiation from hydrogen and cesium 
needs to be included. Another set of parameters needing improvement is the 
viscosity and thermal conductivity in both LHTE and CLHTE. At present, 
the viscosity is that of hydrogen at 30 atn, and the thermal conductivity 


is obtained by use of a constant Prandtl number. The representation of 
these transport properties should be improved and generalized to other 
pressures. Also, the mixing process at the interface between the LSC 
wave heated gas and the buffer gas in CLHTE should be considered more care- 
fully, because of the very large velocity difference. Refinement of the 
stream function grid size should be attempted, to see 'f a non-uniform grid 
IS needed, with closer spacing near the wall, to better define the wall 
boundary layer. The present pressure gradient subroutine provides only a 
constant value for that parameter. Other forms of this function should be 
tried, to learn how to produce desirable channel wall shapes. 



IV. EXPERIMENTAL MEASUREMENTS OF ABSORPTION COEFFICIENT OF 
10.6 ym RADIATION BY WATER VAPOR 

In the continuous absorption uiode, the propellant gas must be seeded 

with molecules which will absorb the laser energy starting at the gas inlet 

temperature, A possible low temperature absorber for 10.6 pm radiation is 

water vapor, which does not dissociate until about 4000 K. It was considered 

as an absorber in the Interim Report (Ref. 4) , and a literature survey was 

made, resulting in a graph of the cibsorption coefficient (in cm ^ per ctnagat, 

19 . 3 

where an amagat is 2,69 x 10 particles per cm ), versus 1000/T. This graph, 
from Fig. 2-5 of Ref. 4 is reproduced here as Fig. 19. Because of the temper- 
ature scale, room temperature is on the right near 3, and temperatures above 
1000 K are on the left, between 1 and 0. 

The dashed line on the left is caused by vibration-rotation line 
absorption and is calculated using a band model, with band parameters based 
on experimental data. The solid lines on the right are caused by continuum 
absorption, and have a pressure dependence, in addition to that from the 
particle density. 

The few measurements of water vapor absorption in the literature up 
to late 1979 are also shown in Fig. 19. There is a group at low temperature 
between 400 and 500 K, and a group between 1535 and 2200 K. The latter points 
are calculated from measured spectral emissivities, and fall nearly on the 
dashed curve. 

Based on this information we have drawn an estimated upper limit 
curve for line absorption, shown by the dot-dash line in Fig. 19. This is 
added to the continuum absorption at low temperature to obtain the absorption 
coefficient used in our modeling. 

Recently, Fowler, Newman and Smith have presented a new calculation, 
and experimental measurements, for the water vapor absorption coefficie.nt 
(Ref. 5). The results are shown as a plot of ku (cm ) vs T in Fig. 20, re- 
produced from Fig. 42 of Ref. 5. The dotted curve is a calculation for a 
specific equilibrium mixture of hydrogen and water, namely = 10/1 at 
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Figure 19 Absorption coefficient of 10.6 Pn 
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Figure CO Calculated and neasurod absorption coefficient of 10.6 .m radiation 
in water vapor fron Fig. 42 of a report by Fowler, Ne\«T:an and Smith 
of United Technologies Research Center (Ref. 5) . 
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11 atm. The values correspond very closely with what would be obtained for 
the water vapor in that mixture by using the dashed line in our Fig. 19. 

The experiments wore performed in a laser heated mixture of 0.3 atm 
water vapor in 13 atm of hydrogen, using a probe laser to measure the absorp- 
tion coefficient. Their results are given by the solid and dashed curves in 
Fig. 20. It can be seen that the measured values presented in Ref. 5 are at 
least an order of magnitude higher than those which they (or we) calculate 
between 1000 and 5000 K. 

This survey of the literature on 10.6 pm absorption by water vapor 
shows two notable features. One is the minimum in the absorption near 500 K, 
which indicates undesirably small absorption coefficients at temperatures 
near the likely inlet gas temperature for a thruster. The second is the 
scarcity of measurements of the absorption, and the scatter near 500 K. 

The only measurement made by actual absorption of a laser line was that of 
Fig. 20 (Ref. 5) , and it shows remarkably high values. 

In view of the scarcity of good experimental results, and of the 
possible importance of water vapor as a molecular eibsorber, PSI has undertaken 
to measure the absorption coefficient in a shock tube. This has the advantage 
of providing a sample of hot water vapor which can be easily prepared, and 
whose properties can be defined with considerable precision (and considerable 
confidence) . The absorption was measured using a probe CO^ laser emitting 
radiation at the P20 line, ilcasurements were made both in the gas behind the 
incident shock and in the gas behind the reflected shock. 

The measurements reported here arc of a preliminary nature, represent- 
ing the first data obtained from the expcrir-ent. The details of the experi- 
mental work are prosei '■od in Appendix I. 

The results are shown in Fig. 21, where the absorption coefficient in 
cm ^ ama ^ is plotted against T, rather than 1000/T, to spread out the region 
of interest. There are 10 points presented, representing the incident and 
reflected reasure-ents on five shots in mixtures of argon and water vapor. 

The arrows on the four left-most points indicate that the lower limit of the 
error bars is not known; indeed, the three left-most measurements yielded 
upper bounds only. The solid curve is a reproduction of the dashed curve 
of Fig. 10, the values currently being used by PSI in modeling. The 
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Figure 21 Measurements of the absorption coefficient of 

10.6 pm (taken at the P(20) line) in water vapor 
from the present experiments. 
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dashed curve of Fig. 21 is drawn at 60% of the solid curve. It falls approx- 
imately through the mean of the data. 

Although the preliminary nature of these measurements must be em- 
phasized, as discussed more fully in Appendix I, indications are that the 
experimental results fall close to the calculations taken from the literature 
in the range 1100 to 3100 K. There is no evidence of the large absorption 
coefficients reported in Ref. 5, as shown in Fig. 20. 

No change has been made in the values of water vapor absorption 
coefficient used in our modeling as a result of these experiments for several 
reasons. First, the results became available only at the end of the present 
effort. Second, the relatively small difference between the experimental 
results and the values being used, combined with the preliminary nature of 
the experiments, did not seem to warrant any changes as yet. 

The experiments are being continued under other sponsorship, and 
when the work is completed, a more informed decision can bo made as to the 
need for changing the absorption coefficients currently used in our computer 
models . 
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V. CONCLUSIONS 


The research study, whose results are presented in the Interim 
Report (Ref. 4) and in the present report, had as its goal the development 
of computer programs useful for the design of thrust chambers for CW 
laser-heated thrusters. To accomplish that goal we have proceeded in 
three steps. /' 


First, we have described the physical properties of some gases ^^ich 
are candidate propellants and seed gases. The properties include laser 
energy absorption properties, broad-band radiation and absorption proper- 
ties, and thermodynamic and transport properties. In connection with the 
absorption of 10.6 pm radiation by water vapor, we have performed a pre- 
limirary shock tube experiment to measure the absorption coefficient. 


Second, we have developed simplified models of thruster chambers 
for boch the continuous mode and the LSC wave mode of laser heating. These 
models can be used for approximate but rapid parametric studies of laser- 
heated thrusters. 


Third, we have developed detailed axisymmetric flow codes useful for 
thrust chamber design for both modes of laser heating ]ust mentioned. These 
codes contain heating by absorption of laser energy, radial energy transport 
by conduction, and radiation energy transport and cibsorption. They have 
been delivered to NASA/Marshall Spaceflight Center and installed on a 
computer there. 

The detailed conclusions which follow from the work accomplished 
are listed below. The first eight conclusions are repeated here from the 
Interim Report (Ref. 4) so that a complete list is collected in one place 
for easy reference. They refer to work reported in the Interim Report. 
Conclusions 9 through 12 refer to work reported in the present report. 
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1 . 


For 10.6 Vtai radiation, a mixture of 1 part cesium 
to 99 parts atomic hydrogen substantially increases 
the absorption coefficient at low temperatures. 

The value of 0.1 cm is reached at 2000 K at 100 atm, 
3000 K at 10 atm, and 12000 K at 1 atm. 

2. For 10.6 Vtn radiation, water vapor has a minimum 
absorption coefficient at about 500 K, and appears 
to be quite low below 1000 K. However, the values 
in the vicinity of 500 K are not well-known, and 
the theoretical predictions have not been checked 
up to 1500 K by experiments. Experimental measure- 
ments between room temperature and 1500 K are needed. 

3. For 5.3 pm radiation, water vapor is a good absorber 
from 300 to 3000 K. NO is a good low temperature 
absorber, and will not dissociate even though it is 
not an equilibrium state. CO is a good absorber 
above 2000 K, and remains chemically stable up to 

at least 6000 K. 

4. For 10.6 Pm radiation, a mixture of up to 5% water 
vapor in hydrogen does not provide good low temper- 
ature absorption. 

5. For 5.3 Pm radiation, both NO and water vapor up to 
5% in hydrogen provide adequate absorption at low 
temperatures, though NO is the better absorber. 

6. The continuously overlapping line model of molecular 
radiation provides results in good agreement with 
other calculations and with experiments. 

7. Cesium is a strong radiator, and its radiation will 
bo an important contributor for thrusters of opera- 
tional size (10 cm radius). It is less important 
for laboratory size thrusters (1 cn radius) , but 
still strong in certain bands. 

8. The thermal conductivity of cesium-seeded hydrogen 
IS nearly the same as for pure hydrogen. The vis- 
cosity, however, increases considerably because of 
the high molecular weight of cesium. 

9. Some limited parametric studies have been made with 
the simplified models of thrust chambers. They show 
that the average wall neating rates are within a 
range that can be dealt with by current rocket 
practice, and that the laser powers needed ar.e not 
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much larger than the power in the exhaust gas. 

At a low power level (30 kw CW 10,6 Ita laser) 
they indicate that small-scale thrusters with 
measurable thrust ran be tested. 

10. The thrust chc.r’ocr design programs produce results 
which are physii.’lly reasonable. However, they 
have only been i ea for a few sample calculations, 
and further expej lep'-e with them is needed to 
oxploio their us^ . 'ness. 

11. Several areas of the modeling in the thrust chamber 
design programs need improvement. The radiation 
modeling for high temperature hydrogen-cesium mix- 
tures does not include cesium radiation at present, 
and cesium may be an important radiator ever in 
small amounts of a percent or so. The representa- 
tion of the transport properties of the seeded 
hydrogen mixtures needs to be improved. The mixing 
process between the LSC wave heated gas and the 
buffer gas should be modeled more accurately. 
Grid-refinement studies should be made to see if a 
non-'.''iform grid in the radial direction is needed. 
The effect of the i"put oressure distribution on 
the resulting channel wall shape should bo studied. 

12. The preliminary experurents on the absorption co- 
efficient of water \apor for 10.6 fim radiation in- 
dicate values which are about 601 of the values 
used in the computer proarams, in the temperature 
range from 1100 to 3100 K. This is considered at 
present to be good confirmation of the values used, 
which were obtained from the literature. The ex- 
periments are being continued at PSI under other 
sponsorship. As yet, no change in the modeling 
has been made based 'n the experimental results. 
When the ongoing experiments are completed, t."'e 
results should be re' lewed to see if -he values 

in the computer prog’-i.ms neoa changinc. 
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APPENDIX A 


MODEL FOR CONTINUOUS HEATING MODE MAPS 

It IS important to find a simple model of the rocket performance 
and wall heating characteristics of a laser-heated thruster. Such a 
model will enable us to quickly survey the range of operating parameters 
and define in an approximate way the combination of parameters which 
may lead to useful thrusters. 

Among the parameters of interest are the thrust F, the- specific 
impulse, Igp» the size, the pressure level, and the average wall heat 
flux. 

In this appendix we will consider the continuous heating mode, 
in which the laser is absorbed continuously over the length of the heating 
zone. We model the heating zone as a cylinder of length and radius 
R^, through which a gas at uniform temperature Tj^ and pressure p is 
flowing at speed u. This cylinder represents the heating zone in a 
laser-heated thruster, as shown in Fig. A-1. We take the temperature 
to be the value after the laser has been absorbed which means the wall 
heating will be over-estimated, so the results will be on the conserva- 
tive side. We tcike the gas to be in thermochenical equilibrium, so 
the pressure and temperature define its density p and enthalpy h, once 
the gas mixture is defined. 

The mass flow rate is 
-> 

m = p u T (A_l) 

the total enthalpy is 


H = h + u^/2 


and the flowing gas is 


P^ = mH 




3U3 


(-< L ^ >.| 



• Conserve mass & enemy 

• Momentum conservation -*■ o constant 

• fisn ■ 

• L = 1/I<L 

• OpTOT ^ ^ 

• Q„ = (Qp ^ Qc)/27tRL 


Figure A-1 Simplified Continuous Heating Model 
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The thrust that can be extracted from this gas, if it is expanded 
to zero enthalpy i . 


F = m u (A-4) 

ex 

where the exhaust velocity at zero enthalpy is 

u^^ = (A-5) 

The corresponding specific impulse is 


I 

sp 


u /g 
ex 


(A-5) 


The hot gas will emit radiation at a rate part of which is 

to the side walls, say Q . It will also transfer heat to the side walls 

H 

of the cylinder at a rate Q^, so the total heat rate to the side walls 

IS 


Q ,1 = Co + (A-7) 

*wall R C 

and the average wall heat flux is 

®w " ^ 2 jJ/(2tt L^) ■ (A-8) 


To maintain a steady state, power must be supplied by the laser, to 

provide the gas power P , the total radiated nower Qp and the convective 

G * ”tot 

power Q_. Therefore, the laser power P, is 
C U 


= 


P +0 
G '^R 




tot 


(A-9) 


The gas power depends only on the thrust and specific impulse, as can 
be seen from Eqs. (A-3) , (A-5) and (A-4): 


D 

*G 


9 


I F/2 
sp 


(A-10) 
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If I and F are given, P_ is not dependent on the size of the thruster, 
sp G 

The energy losses, on the other hand, are dependent on size. The 

convective heating can be estimated from laminar boundary layer theory, 

since these are fairly low-speed flows with moderate Reynolds numbers. 

The formulation of compressible boundary layer theory is in terms of 

the gas-stream density p and viscosity u» end the density, viscosity, 

thermal conductivity, and specific heat at the wall, p , y , X and c 

WWW pw 

The variables used are Q cind n, which are related to the distances along 
and normal to the wall, x and y, by 

X y 

C = .^pyudx, r\ = p\i J lp(x,y)/pldy/(2^)^^^ 

0 0 

where p (x,y) is the local density in the boundary layer. 

In these terms the local convective heat flux to the wall is 




X p uH 
w w 


c (2^) 
pw 


1/2 



where g = H(x,y)/H, the local to stream total enthalpy ratio. The gradient 

(3g/3n) must be obtained by integration of the^boundary layer equations, 
w 

For the present case of constant properties in the axial (x) direction, 
the expression for ^ can be integrated so becomes 


^ = 



w_ 

Pr 



(A-11) 


where Pr = c U/^, the Prandtl number. 

P 

The quantity in square brackets comes from boundary layer theory, 
and docs not depend on x in the present uniform flow. If we let 


G' 


P U 

w w 

p y Pr 



(A-12) 


we can then integrate Eq. 
and lenqth to get 


(A-11) over the cylindrical surface of radius 


"h 
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(A-13) 


= J" ?.TT dx = 2/2? G' H/VJ n) 

0 

This shows that is proportional to i/u, since m is proportional to u and 
H is nearly independent of u in these subsonic, hot flows. Thus, van- 

ishes, i.e., m = 0. IS also proportional to H = h, and increases as 
the flow temperature is increased. 

A suitable value of G* can only be found by referring to boundary 
layer calculations. We are interested here in a gas composed mostly of 
hydrogen, and hydrogen boundary layer calculations have been reported in 
Ref. A-1. Based on these calculations, a good average value of G' is 0,7, 
which means that Eq. (A-13) yields 


= 3.5 H/VJ A 


(A-14) 


This IS the expression we have used in the map calculations performed during 
the present study. 

The radiation emitted by the gas depend both on the geometry of the 
gas Scunple, and on its radiative properties. The radiative properties depend 
on the gas composition, pressure, and temperature, and are expressed in terms 
of an absorption coefficient k^, which depends on the frequency of the radia- 
tion, V. We have studied k^ for gases of interest. For temperatures up to 
6000K, information on k^ is provided in Chapter III of Ref. A-2. For hydro- 
gen at higher temperatures, the work reported there has been extended, and 
is reported in Appendix C of the present report. 


The geometrical dependence of the radiation required here is associ- 
ated with a uniform gas cyli.nder of radius R^^ and length 

total radiation emitted, , and that portion emitted out the cylindri- 
cal surface. Appendix D is devoted to an approximate calculation of these 
quantities in terms of k^. The results are expressed using k^^, R^^, and 
the black body function 


hi- 


We need tre 


B 



exp(h- Ci/kT)-l 

P 


(A-15) 
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where h^, k and c are the Planck constant, the Boltzmann constant, and the 
speed of sound. Rather than use the frequency V, we use the wave number 
0) •» v/c. The result for total radiation is 



For radiation to the side walls, the expression depends on the cylinder aspect 
ratio. 




5 


It 



(1 . VH,’ 


(A-17.) 


00 

0 

4 

= 2ttRj^ c (£) o t 

G = O.TSd-K^) [^1 - Injie/d-K^)] 

+ 1 - 2.25 (I-K 2 ) + 0.25 (I-K 2 ) 

The integrations are performed over the a) range for which B is significant. 

5-1 

For T £ 6000K, the range used is u = 1000 to 10 cm 

We can now describe how this model is used for an approximate cal- 
culation of a laser-heated thiu.ter in the continuously heated node: 

1) Specify the gas .mixture, pressure and , temperature of interest. 
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2) Calculate the gas enthalpy, density and chemical composition 
from a thermochemical equilibrium program. (The 1^^ to be found below will 
be approximately /5h. ) 

3) Choose the length of the cylinder. In the calculations made 

here we have chosen Lj^ =• absorption coefficient of the 

gas mixture at the pressure and temperature chosen, for the laser wave-length 
of interest. This is certainly a characteristic length of the absorption 
zone,' but it is the shortest such length, since the temperature chosen is 
that of the fully heated gas, and k increases with T. A possible refine- 
ment to the calculation is to choose some combination of the absorption co- 
efficients of the inlet gas and tne high temperature gas. The values of k 

la 

are found by using the expressions given in Chapter III of Ref. A-2. 

4) Choose a dcsiicd radius and thrust F. 

5) Find the gas flow speed from Eqs. (A-1) , (A-2), (A-4) and (A-5) . 

t 

These provide a quadratic in u“ whose solution is 
2 r 2 ■’T I /2 

u = -h + h + (F/n 71 (A-18) 

6) Find H, m, P , I from Eqs. (A-2), (A-1), (A-3) and (A-6) . 

G Sp 

7) Find the heating terms Q_, Q , Q from Eqs. (A-14) , (A-16) 

^ t ^ 

and (A-17) . In Eq. (A-14), the gas viscosity is approximated as that of 
hydrogen, which is given in Table A-1, as calculated by the method described 
in Chapter IV of Ref. A-2. As pointed out there, small fractures of heavy 
seeds like cesium can increase the viscosity of the mixture by as much as 
a factor of 2 over pute hydrogen. However, the convective heat transfer 
term is already approximate because of the choice of G’ , and it only depends 
on the square root of j, so for this simple model, it was not felt worth 
while to consider the dependence of U on gas composition. 

8) Find the wall heat load average wall heat flux 

0 from Eqs. (A-7) and (A-8) . 

w 

9) Find the laser power P from Eq. (A-9). 

Lf 
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TABLE A-1 


VISCOSITY OF EQUILIBRILM HYDROGEN 
(in 10 ^ kg/m-s) 


T(IO^K) 

Pressure (atm) 

1 

3 

10 

30 

100 

1 

1.98 

1.98 

1.98 

1.98 

1.98 

2 

3.13 

3.13 

3.13 

3.13 

3.13 

3 

4.48 

4.42 

4.39 

4.37 

4.36 

4 

5.27 

5.58 

5.68 

5.64 

5.55 

5 

5.54 

5.78 

6.25 

6.77 

7.13 

6 

6.54 

6.59 

6.74 

7.06 

7.62 

7 

7.33 

7.35 

7.40 

7.55 

7.94 

8 

8.32 

8.34 

8.38 

8.45 

8.66 

9 

9.10 

9.18 

9.23 

9.30 

9.42 

10 

9.72 

9.95 

10.1 

10.2 

10.3 

11 

9.97 

10.5 

10.9 

11.1 

11.2 

12 

9.54 

10.6 

11.3 

11.7 

12.0 

13 

8.49 

10.2 

11.5 

12.2 

12.7 

14 

6.95 

9.30 

11.3 

12.5 

13.3 

15 

5.25 

7.93 

10.5 

12.3 

13.6 

16 

3.76 

6.49 

9.66 

12.1 

14.0 

17 

2.56 

4.91 

8.21 

11.1 

13.6 

18 

i.ao 

3.63 

6.74 

9.96 

13.0 

19 

1.32 

2.70 

5.44 

8.75 

12.3 

20 

1.04 

2.06 

4.34 

7.54 

11.4 

21 

0.881 

1.64 

3.48 

6.41 ■ 

10.4 

22 

0.824 

1.42 

2.85 

5.44 

9.43 

23 

0.776 

1.26 

2.40 

4.67 

8.52 

24 

0.756 

1.15 

2.08 

4.05 

7.69 

25 

0.754 

1.10 

1.98 

3.58 

6.97 
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This procedure has been implemented in a smaxl computer program which 
uses the equilibrium state and k as input. Calculations have been made for 

U 

a number of cases. The results are presented in the main body of this report. 
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APPENDIX B 


MODEL FCA LSC WAVE HEATING MODE- MAPS 

This Appendix describes a simple models of a laser-heated thruster 
which uses an LSC wave mechanism to absorb most of the energy of the laser. 
Since that mechanism produces a-high temperature-gas— (-20^000K_in^-pure 
hydrogen), it is hard to contain by walls, and produces a great deal of 
radiation. To make it useful as a thruster, we consider surrounding the 
heated gas by a co-flowing buffer gas of the same mixture, but unheated 
by the laser. This buffer gas will be heated by the laser-heated gas by 
both convection and radiation. The two gas streams will mix and ultimately 
produce a single stream at some temperature between that of the hot gas 
exiting the LSC wave and that of the incoming buffer gas. The buffer gas 
will also serve to produce the convective and radiative heating to the 
walls. 

A schematic diagram of such a heater is shown in Fig. B-1. The 

over-all size is radius R and length L. + L . Flow at state ( ) enters 

the heater. In the center, over a radius a laser beam of power P^ 

enters, inducing an LSC wave in the cylinder of length Lj^. This hot 

cylinder radiates energy some of which escapes out the side to the 

walls, and some of which is absorbed in the cooler buffer gas in the 

annulus R^-R^^. There is also convective heating from the hot gas to the 

buffer gas along the edge of the hot cylinder. Starting at Lj^, the 

gases mix, and at the exit of the heating region they emerge as one 

stream with properties ( In the mixing region, the gases lose 

heat to the side walls by convection, Q , and radiation, Q . Heat is 

On Fm 

also radiated out the ends of the mixing - region, AQ . No radiation out 

Kltl 

the ends of the LSC wave cylinder is considered, as it is all absorbed in 
the oncoming gas and the mixed gas. 

This model permits a straightforward calculation of the flow param- 
eters and losses to be made by using conservation of mass, momentum and 
energy. Actually, momentum conservation is replaced by the condition of 
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» VL3 




^Rh “ ^wall . 

• Conserve mass s eneray 

• Momentum conservation -vp constant 

r 9 1 1/2 

• °^sp “ Pmf ^ %] 

• ^Ti ■ 2 Rjp^ L|^ = l/*S/iiv 

• ”1 ■ 8Rr.tot * 

• ^ - (Qpi^ + >• Q^p^)/2 7t R j^L^ 

t I[_ related to p and m^ 

Figure B-1 Siinplified LSC Wave Keating Model 
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constant pressure throughout, leaving the mass and energy equations to pro- 
vide us with the relations necessary for the calculation. 

Let us denote the gas conditions at the end of the LSC wave region 
(hot region) by ( and those at the end of the buffer regie ,i by ( )^. Then 
over-all mass conservation requires that 


ii + ITL = p. u S 
b n 1 1 n 


p u T 1 R 
11 m 


P , u . s 
mf mf m 


(B-1) 


where is the cross-sectioftal area of the cylinder. The stagnation enthalpy 
at the end of cylinder is 


H = h + u /2 
mf mf mf 


(B-2) 


and the power in the gas is, with the input gas power, 


P^ = mH^-P, P =mH 
G mf 1 1 1 


(B-3) 


The thrust resulting from expanding this flow to zero enthalpy is 


F = m u 

ex 


(B-4) 


where the exhaust velocity is 


u 

ex 


(2H ,) 
mf 


1/2 


(B-5) 


and the corresponding specific impulse is 


I = u /g (B-6) 

sp ex 


These relations are the same as those in Eqs. (A-l)-(A-6) for the 

thruster performance, once the gas is heated and mixed to the ( ) , state. 

mr 


The radiation associated with the mixing region flow is composed 


of the total emission o 


Rm 


, the radiation to the cylinder side wall 


tot 
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■» Qj^ - These are calculated from the pressure, temperature 

and cylinder geometry R^, by the formulas given in Appendix A, 
Eqs. (A-16) and (A-17), using the appropriate gas composition. 

The convective heating is found from Eq. (A-14) of Appendix A, 

using H L , and the viscosity u , at T ,, according to 
mf m mf mf 


2cm='-' 


H , /p , m 
mf mf 


(B-7) 


To obtain information about the conditions at this exit of the hot 
region we write an energy balance between the station and the end of the 
cylinder. The energy flow emerging from the hot regio.i is m^^ That 

emerging from the buffer region is m^^ Here is the 

convection loss from the hot region, which we take to be deposited n ’■he 
buffer gas. It is given by the formula analogous to Eq. (B-7) a”" 


«ch- 

Q IS the energy lost to the wall by convection from the buffer gas, which 
Co 

we will ignore as small, since the buffer gas is much cooler than the hot 

* 

region gas. Finally, is that portion of the radiation from the hot 
region which is absorbed in the cold region. It is found by subtracting 

from the radiation emitted by the hot region, Q„. , that portion which can 

■ 

escape through the buffer gas to the walls, Qpj^* 

For hydrogen, these quantities are found as follows. Q , comes 

Ril 

from using Eqs. (A-17) , with the integration on m carried as far as necessary 
to obtain significant contributions; in all calculations we have integrated 
to t) = 2 X 10^ cm with k^ taken from the high temperature hydrogen radia- 
tion properties described in Appendix C. The part which can escape through 

t’->'“ buffer gas is calculated the sane way, but only integrated up to 
4 -1 

= 9 X 10 cm , since the higbei energy part of the radiation is absorbed 
i.n the cooler buffer gas. 

The sum of the energy flows emerging from the buffer and hot regions 
either emerges from the end of the mixing region as m or- goes to 
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the walls as radiation from the mixing region, Q , or convection Q 

Rm. . Cm 

tot 

The energy balance from station to station Lj^ + is then 


m.H,+m,H +Q^ + Q_. = m H , + Q_ + Q 
n li 1 *'ch Rh mf ^ 

tot 


(B-9) 


This can be considered a quadratic equation for since m^ = m - lii^, 

and Q depends on , Eq. (B-8) . Its solution is 

Ch n 


. 1/2 


2(Hh - H^) 






* 1/2 

m(H , - H ) + + Q - Q ^ \ 

' mf i' f^tot I 

Ku - H j 

n 1 ' 


(B-10) 


We can now set up a procedure for calculating the performance of an 
LSC wave heated thruster with a buffer gas. The parameters to be specified 
are partly determined by convenience of calculation and partly by the prop- 
erties of LSC waves. It is known (Ref. B-1, for example) that there is a 
i elation between laser intensity ard mass flow rate per unit area for LSC 
waves. This relation should be maintained for the calculations with this 
simplified model. Therefore, we choose the input mass flow per unit area 
u^ as a parameter specified in the calculation. In order to make 
the thermochemical calculations as simple as possible, we specify p and T 
to determine the state of the gas, wherever we need it. With these con- 
straints, the calculation proceeds as follows: 

1) Specify the gas mixture, pressure p, thrust F, inlet tempera- 
ture T , mass flow per unit area .m. , and final mixed temperature T 

1 A mf 

2) Calculate the enthalpy h , density p ,, and chcnical comoo- 

inf mf 

sition from a thermochcmical equilibrium program. 

3) Calculate u . = m,/c ,, H , from Eq. (B-2) , n from Eqs. (B-4) 

mt A mf mf 

and (B-5) , R from Cq. (B-1), I from Eqs. (B-5) '-’d (B-6) . 

m sp 

4) Calculate h^ and from p, using perfect gas properties; 

then u = m./P , and H = h , + u^/2, and the i.iput power P. = m H •. 
lAl 111 r-i- ill 
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5) Choose Ii , the length of the mixing region. This choice can 

m 

be based on a length of mixing, or some other suitable interim. For the 

calculations made in this study we have chosen L = 2 R . 

mm 

6) Calculate the heating terms Q , and from Eqs. 

Om f^n ^uTi 

(B-7) , (A-16) and (A-17) . In Eq. (B-7), the viscosity b , is approximated 

mi 

by the value for hydrogen from Table A-1. 

7) Calculate the power in the gas at the end of the mixing region 
by Eq. (D-3) . 

8) Choose the temperature at the end of the LSC wave (hot 
region) and calcruate the corresponding gas composition from a therroochem- 
ical equilibrium program, including enthalpy h^^ and density pj^. 

9) Find the velocity Uj^ - m 

10) Choose the length of the hot region, This should be a 

length characteristic of the thic)cnoss of the LSC wave. It is not the 

absorption length at Tj^, since this is usually very small. In LSC waves, 

the thic)cnesb is determined by the distance over which the cold incoming gas 
IS heated up to a temperature at which it can absorb the laser energy. This 
heating is caused by conduction and radiation forward from the hot absorbing 

part of the wave. In hydrogen, the length over which this heating occurs is 

characterized by the absorption length for vacuum ultra-violet radiation 
emitted by the hot hydrogen. According to Appendix D of Ref. B-1, this 
length is, in cgs units, 

— 18 

K = 6.34 X 10 (2 r/Q . ) (B-11) 

vuv H ^e«.H 

where n is the hydrogen number density, and Q , is the electronic partition 
H Q tH 

function for hydrogen atoms For the puce hydrogen calculations done in the 
present study, wo have taken Lj^ = 

11) Guess a radius of the hot region. Then calculate the 
radiation from the hot region, as described above, by using Eqs. (A-17) 
integrated over the whole significant wavenumber regime (0 to 200,000 cm ^ 
in the cases calculated here) . This gives the total sideways emission 


and total enthalpy = hj^ + 
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Q . . Then use the s^une equations to find that part of the sideways emission > 
Rri 

Q* , which can escape through the buffer gas to the walls, by integrating 
over the appropriate portion of the wavenumber spectrum (0 to 90,000 cm 
in^the. cases -calculated here). Subtraction then yields »• 
the power absorbed by the buffer gas. 

12) Calculate the mass flow rate mj^ through -the hot region from 
Eq. (B-10) , which then gives the radius of the hot region from iiR^ => 

*> If this radius is not the same as the one guessed in Step 11 

<ibove. Iterate between Steps 11 and 12 until the radius calculated from 
Eq. (B-10) is the same as the one used to find the radiation in Step 11. 

13) Calculate the laser power as the power in the gas at the end ^ 
of the mixing zone, plus the power put into the walls; 


P=P+0 +0 +0, 

L G '^Rm ^ ^Rh 

tot 


(B-12) 


The laser intensity is then 

14) The total wall heat load is 


^all ^Cm ^Rm * ^Rh 


(B-13) 


and the average wall heat flux is 


? = 0 ,,/2n R (U + L ) 

w wall tn h m 


(B-14) 


15) To satisfy the constraint that the input value of is appro- 
priate for a given laser intensity I . we may iterate on any other of the 

Ij , 

* 

input parameters such as Only one need be changed, as only 

one output parameter, I , is being sought. Wo have found it convenient to 

Ij 

Iterate on T. , since the thrust F often has a desired Value, and T , almost 
n tnr 

comoletely determines I , which also often has a desired value. 

sp 
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The procedure through step 14 has been implemented in a small com- 
puter program for pure hydrogen propellant. It takes as input the equili- 
brium properties of hydrogen at the two states p, and p, Tj^. The iter- 
ation on I , or any other parameter, is done manually using multiple runs 
of this computer program. 

A number of calculations of typical cases have been performed with 
this program, and the results are discussed in the main body of this report. 
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APPENDIX C 


IlYDRfjGEN RADIATION AT HIGH TEMPRATURE 

Radiation Cron hydrogen species is discussed in Chapter III of Ref. 

C-1. However, the expressions given there are limited to a temperature 
regime up to 10,000 K. Tor the case of LSC waves, especially in pure hydro- 
gen, tenporaturos up to 20,000 K, or somewhat higher, are of interest. Such 
waves are modeled in the hot regime of the simple model described in Appen- 
dix B, as well as in the detailed structure of hydrogen LSC waves described 
in Ret. C-2. Therefore, wo need to extend the radiative properties of hydro- 
gen to these higher temperatures. 

There are two effects which enter as higher hydrogen ten^eratures 
are considered. First, the significant part of the frequency spectrum ex- 
tends to higher values, which is to say, to higher wavenumbers, because of 
the shape of the black body spectrum. As temperature increases, the peak of 
the blackbody curve shifts to higher frequency, and its fall-off also shifts 
to higher frequency. For a temperature of 10,000 K, the black body curve has 
fallen to about 10 ^ of its peak value when the wavenumber is 10^. The same 
degree of fall-off at 20,000 K requires a wavenumber of 2 x 10^. Integrations 
over wavenunber, which' were carried out to wavenumber 'j = 10^ in Ref. C-1, were 
extended to i.,' - 2 x 10^ for high temperature hydrogen. 

The second effect of high temperature is the inclusion of hydrogen 
line ladiation, to be added to the continuum radiation considered at low 
temperature. Most lines can be treated as extensions of the continuum radia- 
tion for a simple model, .owever, the Lyman alpha (Ly^) and Dalmer alpha (H^) 
ii.nes may be strong enough to be treated separately. 

The model wo have adopted for the absorption coefficient of high tem- 
perature hydrogen is described oelow. It consists of modifying the free-bound 
hydrogen absorption coefficient given in Ref. C-1 to allow for wavenumbers up 
to 2 X 10^, of adding terms to model the lines as an extension of the continuum, 
and of adding the Ly ^ and lines. 
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We start with the expression for the free-bound absorption coeffi- 
cient of hydrogen atoms given in Eg. (J-20) of Hef. C-1. 


k 


H 


'h "h 


-hpCO)/kT 

Q 


) 


(C-l) 


Here 0„ is the cross-section, n„ the atom nuiiijer density, h the Pl£mck con- 
H H p 

stant, c the speed of light, w the wavenumber and k the Boltzmann constant. 
The cross-section is expressed in the equation above Eq. (3-20) in Ref. 

C-l in cgs units, as 


o„ = 


2 X 7.9 X 10 
^ciH 


-18 


max 

E 

N > /Ry/u) 


3 

/Ry\ 

r •'y / 1 V 

(m) 

kT r 2 ) 


L \ N /J 


where Ry = 109737 cm ^ and the electronic partition function of the hydrogen 
atom IS 


N 

m.ix 
N=1 

IS the largest principal quantum number of an electronic state of hydro- 
gen to be used in the sums. At low tc-’peratures, the choice of this cut- 
off quantum number does not ratter much, since the larger N terms make little 
contribution to the sums. Hoi«ever, as T increases, this ceases to be true 
and a realistic cut-off must be chosen. While in the low temperature work 

of Ref. C-l we chose N “ we have chosen N = 10 in the present 

max max ^ 

high temperature case. A more sophisticated cut-off procedure was not con- 
sidered worthwhile for the present rat.~.er sir^le model. 

The first modification to Eq. (C-2) is for the N = 1 term, which 
never arose when < 10^ < Ry, since then the lower limit was always N > 2. 


hpC Ky 

/l _ 

kT 

V nV 


(C-3) 



But if we extend the range of w to 2 x 10^, 

the higher wavenuznbers . This term is given 

the last term on the right. In the present 

ground state number density n,_ = 2n„/Q o.,* 

-18 " S/3 

sum in Eq. (C-2) is 6.34 x 10 {Ry/<jJ)° 


the N = 1 terra will appear for 
i.i Ref. C-3, Eq. (5-17), as 
units 13. 6/(1) = Ry/U), and the 
Thus the N = 1 term in the 


as 


To add the linos as an extension of the continuum we write Eq. (C-2) 


10 



N=1 


(C-4a) 


The A , are chosen at each (d as follows. Define 
N 

= Ryd - 1/9) = 8 Ry/9 

OJ, = Ry(l/4 - 1/16) = 3 Ry/16 (C-4b) 

N>2: = Ry (n”^ - (N + 1)“^J 

Then the A are; 

N 


N = 1: 


0 if (i) £ 0)^ 

6.34x 10 (Ry/^)^^^ if u '• 'i), 


(C-4c) 



To conpleto the model, we need a treatment of the Ly^ and lines. 
Those lines arise from the transitions n = 1 to n ' =2 and n = 2 to n’ = 3, 
respectively. The line centers are at given by 


‘'l2 = = 3Ry/4, = Ry.il/4 - 1/9) = 5Ry/36 



which, we may note, are values missing from the sequence in Eq. (C-4b) , 
which have the n = 1, n' =• 3 and n = 2, n’ » 4 transitions. The correspond- 
ing f , numbers are 
^ nn' 


12 


0.4162, 


^23 


0.641 


We first find the absorption coefficient integrated ovei the lino 


width as 


fk d(i,= 

J to 




4 exp 


/ A 

\ kT 4/ 


where the upper and lower factors in { } are for Ly^ and H^, respectively. 

Here e is the electronic charge, m the electron mass and n is the number 

£< ACi 

density of atoms in the ground state. The latter is the lower state for 

the Ly line. The lower state for the H line is the n = 2 state, which 

™ 2 “2 
requires for its number density the added factor 2 exp [- h^c Ry/kT(l - 2 )] 

which explains the factor in { } for H^. In terms of the hydrogen atom 
number density we have n^^ = 2n^/Q^.^, where Qqjjjj the electronic parti- 
tion function for H, already given in Eq. (C-3) . Using this expression for 

n.^ and evaluating the constant, we find in cgs units 
AG 


fk d.' = 


2 X 8.S5 X 10“^^ f , n„ 
nn H 


Q 1 


I - exp 


I V 

\ ) 


4 exp 


I 3 

\ 4 kT / 


(C-5i 


We now evaluate the continuum absorption coefficient at the line center, 

say k ^ (CJ , ) , by the formulas given for hydrogen in Chapter III of Ref. C-1. 

^ cont nn' 

Let the gas sample length we are interested in be L. If k ._(:j ,)L > 1, 

’ cont nn' 

the continuum contribution to emissivity is already near the black body limit 
and the addition of the line contribution to the absorption coefficient does 
not increase the emissivity significantly; therefore the line contribution 
can be neglected. 



If k ^(u) ,)L < 1, we then proceed to find the line contribution, 

cont nn ' 

One possible value of emissivity for the line is 


«= B(o) ,) L f k doj/oT'' 
nn' J (jJ 


(C-6) 


where a is the Stofan-Boltzmann constant, B(u) is the black body function 


B(w) = 


^ V 2 3 

27ih_c w 
P 


exp(hpC (jj/kT) - 1 


and the integral of the absorption coefficient is given in Eq. (C-5) . This 
value IS appropriate if the line is transparent, i.e., if the emissivity at 
the line center is not limited by reabsorption. The transparent approxima- 
tion is used whenever [k (oj .) + k ^(o) ,)]L < 1. 

(i) nn* cont nn' 

If the line center is reabsorbed, the line emicsivity can be approxi- 
mated by using the transparent emissivity in the spectral region where 

< 1, and the black body limit in the spectral region about the 

line center where (k ^ (u) ,) + J: ]L > 1. This second possible emissivity 

cont nn' w ^ ' 

IS obtained using a line width 


Au' 


B 


^ 2 
2a . 

T nn* 


F 

o 


where (in cgs and vtatcoulombs) 


F = (4Tjn = 1.248 x lo”^ 

o E E 


and a„ is found from 
T 


C , /k d. 
nn ' ■' 


. F 

nn o 


Fl"^ - k ^ U< , )]\ 

L cont nn' Jl 


2/d 



I. 
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The constants C . for the two lines are 
nn' 


= 3-4 X 


C 23 “ 1.3 X 10 


-15 


in consistent xinits. Then the second cmissivity is 

.( 1 ) 


.( 2 ) 


D (oj , ) o 

nn' B . 8 


'»T> 


3/2 


(C-7) 


where the first tern represents the contribution of the black "line center and 
the se"ond term is the contribution from the treinsparent line wings. 

The emissivity contribution of each line (Ly^, n = 1, n' = 2 and 
n = 2, n’ = 3) is given by the smaller of and for the line. 

This IS added to the continuum contribution to e, which is 


e 


cont 



(1 - e"’'^) du 


(C- 8 ) 


where k is the sum of the contributions discussed in Ref. C-J., Chapter III 
for hydrogen species, but the free-bound absorption coefficient of hydro- 
gen atoms given there is to be replaced by that given in this Appendix, 

Eqs. ;c-l) - (C-4). To avoid counting the continuum twice in the spectral 
region where the I’jie center is rccibsorbed, we subtract from the C for the 

lines the value of the integrand of (C- 8 ) at (o , times A£i)„. 

nn' B 

The cmissivitics that result from this model are displayed in Table 
4 4 4 

C-1 for T = 10 , 2 X 10 and 3 x 10 K, p = 1, 10 and 100 atm, and path 
lengths of 1, 10 and 100 cm. The inc’-eases with pressure and path length 
at a fixed temperature is shown clearly. The values cover a large range, 
from 10 ^ to 1 . 

A previous calculation of the emissivity for high temperature hydro- 

4 4 

gen has been published in Ref. C-4, for 10 and 2 x 10 K, with 1 cm at 1 
and 10 atm, and 10 cm at 1 atn. The values of emissivity given there are 
about a factor of 2 higher than the ones we have calculated. However, our 
description of the line shapes is more modern than that of Ref. 0^-4, which 
probably accounts for the differences. ' » . 
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EMISSIVITY 


TABLE C-1 

HIGH TEMPERATURE HYDROGEN 



2. = 10 cm 

1 = 100 cm 

1 atm 

- 


-3 

-2 


5.57 X 10 

3.74 X 10 


-2 

-1 


5.67 X 10 

1.84 X 10 


-3 

-2 


5.63 X 10 

4.62 X 10 

10 

atm 



-2 

-1 


6.03 X 10 

4.12 X 10 


-1 

-1 


5.00 X 10 

9.94 X 10 


“1 

-1 


2.80 X 10 

6.46 X 10 

100 atm 


-1 



6.75 X 10 

1.03 


-1 



9.86 X 10 

1.00 


-1 

-1 


9.32 X 10 

9.86 X 10 
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APPENDIX D 

RADIATION' FROM A CYLINDER 


The simple thruster models described in Appendices A and B require 
an approximate expression for the radiation from a cylinder of gas at uniform 
temperature and pressure. The purpose of this Appendix is to derive the 
expression used in those models. 


Consider a surface enclosing a uniform body of gas. The radiative 
flux normal to the surface at any point of the surface is made up of the 
contributions of all the rays which can reach that point from within the gas. 
Let each such ray bo identified by angular coordinates 0, (j>, where 0 is the 
angle the ray makes with the inward normal, and ip is tlie azimuth angle around 
the normal. Let be the length of the ray, from the surface point to 

the place where it pierces the surface again on the other side of the body. 
The radiative flux normal to the surface, at frequency v, (and wave number 
o^v/c) is expressed by an integral over all such rays, as 


I f 

U V J 


2tt 

0 


tt/2 


d0 sinO COS0 


1 - z 


.k^i(6,4.) 


(D-1) 


where k 

(J 


IS tlie absorption coefficient at o). 


is the black body spectral flux 


D 

U) 


2 3 

27Thp c w 


exp (hpCu/kT) - 1 


(D-2) 


and hp, k are the Planck and Doltzmam constants respectively. 

To find the power radiated from the whole body at wave number u) 
we must integrate over the whole surface 

= / S dA (D-3) 

\ot 

and similarly, to find the power radiated from any portion of the surface 
we perform the integral over that portion. 



Finally, to find the total povcr, wo must integrate over U 



dOJ O 

’^toto) 


(D-4) 


and similarly for any portion of the surface. 

It 13 thus apparent that for a unifon body of gas, the calculation 
of the radiation can be thought of in two parts. One is the gcome‘-ric prob- 
lem of calculating the integrals (D-1) and (D-3) . The other is the physics 
problem of defining as a function of gas- composition, pressure and tempera- 
ture, and integrating over (j. These two parts can be trcatca independently. 

In the present work the integral over O) is performed numerically, once k^ is 
cetcrmined. The present Appendix is devoted to the geometrical problem. 

The integral in Eq. (D-1) is not expressible in closed form for most 
bodies, but can bo performed for some simple shapes. For a plain infinite 
slab of thickness d, the classical result is 


-= 1 - 2E3(k^d), E3 (x) 
0) 


oo 

r -tx d 

7 ° ; 


d^ 

3 


(D-5) 


For a sphere of radius R wo easily find 


- = 1 + 


-2a 


:(l-e"“^) 

~> 


, a = k R 
0) 


(D-6) 


For an infinite cylinder, the result has recently been giien (Ref. D-1) in 
toins of ncdified Fosscl functions as 


-= ^ |:[aK3(a)l3(a) t aK^(a)I^(a) - ij 

- K^(a)l3(a)| 


(D-7) 


+ K3(a)l3(a)/a + <3(a)I^(a) 




J. i 



All these results are, of course, independent of position on the surface, 
because of symrctry, so surface integration is simple. 

We are here interested in a finite cylinder of radius R and length L, 
for which Eq. (D-1) cannot be integrated in terms of known functions. We will 
therefore consider an approximate expression. 

Tiie limiting forms of Eq. (D-1) for very large and very small kji 
are called the black and transparent limits. The former is obtained by 
letting the exponential vanish, whence we find 


0)D 


= .H 


'jJ 


(D-8) 


as we expect. The transparent limit is the first term of the expansion of 

the exponential, loading to 

, „ 27T 27T 

k Q ^ ^ 

^(OT ’ ' ^ ^ J J (D-9) 

0 0 


The coefficient of is called the geometrical mean beam length 


, ^ 2tt ^2tt 

= — / d^ij d0 sinO COS0 i(0, 
0 


<J» 


(D-10) 


A frequently used approximation for gas radiation is 


S 

U) 


= B 


( 1 ) 



Jo 


) 


(D-11) 


This expression has the correct value in both the black and transparent limits, 
Eqs. (D-8) and (D-9) respectively, and is often used as an approximation over 
the whole range of k 

CJ 0 

Some idea of its accuracy can be obtained by applying this idea to the 
exact expressions for the slab, sphere and cylinder. The values of obtained 
from k d or k R -*• 0 in Eqs. (D-5, D-6 and D-7) , are 2d, 4R'3 and 2R resoectively 

CJ 

These are then substituted into Eq. (D-il) , and the results compared with exact 
evaluations of Eqs. (D-5, D-6 and D-7). (Tables of and the Bessel functions 

are available for these evaluations.) Such a calculation shews that the maximum 

\ 

' t 


' 19 ' 



error for the slab is 13.5% at k d = Q.5, for the sphere it is 5.2% at 

uJ 

k R - 1.5, and for the cylinder it is 6.4% at k R = 0.8. Thus the use of 
hi w 

the approximation fD-11) is quite reasonable for these cases. 

In the same spirit one can intcqrato Eq. (D-10) over the surface 
area of interest to obtain what wo can call an area-averaged geometric moan 
beam length. 

ciA ID-12) 

A 

The limiting radiated powers are then the black value 


^R(i)B 


AB 


U 


(D-13) 


and the transparent value 


Q 


RuT 


k B Ail„ 
(U (1) 0 


(D-14) 


Wo can then propose an approximation for the radiated power, in analogy to 
Eq. (D-11) , as 




= 


RtiiB 




(D-15) 


This again has the correct black and transparent limits, and can be thought of 
as u possible anoroximation over the whole range of k 

Checking the validity of this approximation is not so simple as for 
Eq. (E-II) . V.'c cannot use the slab, sphere and cylinder, since their values 
of do not vary over the surface, so Eq. (D-15) reduces to Eq. (D-11) times 
A check is possible in- a fairly simple way for the case of the base of a 
scmi-inf inite cylinder, for which it can be shown that the exact expression is 


- 1 / -i- 


(D-16a) 
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Tho expression obtained for 2^ is 16R/3n so the approximation is 




Qr = hr B 
01 


(D-16b) 


A niimct leal integration of Eq. (D-16al shows that Eq. (D-16b) has a maximum 
error of 27,5% of k R = 1. This gives us some oonfidonco in using Eq. (D-15) . 

(J 

Certainly other chocks would be desirable, but wo know of no other exact ex- 
pressions which have boon evaluated. Such evaluations would be very helpful. 

It IS felt that present evidence shows Eq, (D-15) to bo sufficiently 
accurate for use in simple models such as those described in Appendices A 
and B. We will therefore proceed to apply Eq. (D-15) to a finite cylinder. 

If one IS i.itorosted in tho radiated jxjwer emitted over tho whole 

surface enclosing tho volume, with area A , it is not necessary to perform 

tot 

any intoerations to obtain tho limiting values. The black limit is, of course. 


? = A B 

"Ral? ^ tot 01 
tot 


(D-17) 


Tho tr.insparont limit is easily found, because there is no absorption, so the 
radiation from oacli voliine element is emitted. In any one direction a volume 
element dV radiates k B dV/TT per steradian, or 4k B dV in all directions. 

iil 01 01 01 

Tlien the total tianspareiit radiated power out of a lolumc V is 


ROJT 


tot 


Ik B V 
u 01 


(D-18) 


This must agree witli tlie integral of Eq. (D-9) over tlie total surface area, 
which moans that 



;v = 


tot 


tot 


(D-19) 


Tliat this IS correct can be checked for the slao, sphere and cv'mdor, wnero 

is 2d, 4R/3 and 2R, independent of location on the surface. Tho corresfxjnaina 
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2 3 2 

areas are 2A, 47TR and 27TRL and volumes arc Ad> 4 ttR /3 and ttR L, where A 

IS the area of one face of the slab. Thus in all three cases i A^ ^ =» 4V, 

0 tot 

in accord with Eq. (D-19) . 

We turn now to the calculation of for the side surface of a finite 
cylinder of radius R and length L, which is needed for the simple nwylols of 
Appendices A and a. Consider first the mid-point of the side. Recall that 
0 is the angle a ray makes with the inner normal from the mid-point, and 4> 
is the azimuth angle around ’•Jiat normal. We choose (|) = 0 in a cross-sectional 
plane. Then a ray from the nid-point of the side in any direction 0, <t>, 
which hits "the side of the cylinder has length 

Z.^ = 2R cosO/Ccos^G + sin^O cos^cj)) (D-20) 

A ray which hits the end of the cylinder, where <p is between 0 and n, has 
length 


>. = L/2 sinG sin<J 

o 


(D-21) 


Because of symmetry we need only consider <J> in the range of 0 £ (J) £ TT/2, and 
then quadruple the rcsulc. Cf course, 0 only vanes between 0 £ 0 £ tt/2. 

The dividing line between whore we ha\e to use Z and where we must use Z 
IS obtained by setting Z^ = Z^. This is the equation of the circular perimeter 
of the base in 0, <p: 


2 2 2 

2R cost? sinG sin6, = (cos 0 + sin 0 cos‘"<J),) L/2. 

d d 


(D-22) 


We can now express from Ec. (D-10) as 


^0 = f / 

V 


deb 


0 


/■ 


dO sin9 cosG Z 


dO sinO COS0 Z 


(D-22) 


where <|>,(0) is defined bv Eq. (D-22). 
d 
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In those variables tht integration is very difficult because of the 
complicated dependence of and on and 6, as well as the complicated 
form of <}'j{0) • However, there is a set of angles which makes these expressions 
much simpler 

Let tt/2 - ct bo the angle between the ray 0, (J and the generator of 
the cylinder through the mid-point. Project the ray on the cross-sectional 
plane, giving a projection of length i cos a. Then introduce the angle between 
this projection and the inner normal from the side mid-point as 3. These 
angles are related to G, <(i by 

cosQ = cosa cosS , sina = sin<f sinO (D-24) 

and £. , a take the simpler forms 
so 

i = 2R cosS/cosa , 2. = L/2sina (D-25) 

s e 

The relation between G and on the perimeter of the cylinder end, Eq. (D-22) 
becomes 


tana, = L/4R cosS 
d 


(D-26) 


Finally, to transform the integral (D-23) we need the Jacobian determinant 

D0/3a 30/33 


J = 


3<)/3a 3<{>/3S 


to express dGdi{> = |j| dadS. The algebra is lengthy but straightforward, and 
we find 

2 

sinO cosS d0d$ = cos a cosB dndS 

The limits on 3 are from 0 to tt/ 2. The limits on a are from a, to tr/2 for 

d 

the rays which hit the end, and 0 to for the rays which hit the side, where 

a, (S) IS given by Eq. (D-26). 
a 
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t 


The expression for Eq, (D-23) , then becomes 




. SSy"'‘d8 coss/' 

0 0 

2L r r 

^ / dg cosg/ 


cosa 


TT/2 2^ 

jn of J, COS a 

dg cosg/ da 

' cina 


a. 


The first integral on a xs easy. The second can be broken into two parts by 
2 2 

using cos a = 1 - sin a, and the second part can also be integrated. The re- 
sult IS 


dg cosgj^4R cosg sincij “ L cosa^ 


0 

2L r 

+ — / dg cosg / 


tt/2 ^ti/2 

dg cosg / da/sina 


But the term in square brackets vanishes because of the definition' of a. from 

a 

Eq. (D-2G) , leaving only the last double integral. This integral can be inte- 
grated by parts on g to give 


'o-f 


TT/2 


-/’"iff. 


j i X tt/ 2 „ da, 
+ / dg — 

sina I J sina dg 


0 0 


The integrated part vanishes because a^ = t/ 2 when 6 = ti/ 2, as Eq. (D-26) 
shows. Also, from Eq. {D-26), wo find that 

1/2 


4R 


sina, L 
d 


(:fe) 


dg “ 4R 


sing 


(L/4R)^ + cos^g 
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Finally, then, we have 




tt/2 

1 ^ ik f ae 7 ^ 

0 ^ / r„ . __2j 


sin^S 


[ 


(P-27) 


(L/rR) + cos P, 


This can be expressed in terrts of complete elliptic integrals by writing the 
integrand as 


1 1 - 1 * ^ ^ f/L.f 


(D-28) 


Thus we have 


^ 

0 Tib, 


2L 

Tib, 


jr/2 

/ dB 


(l - b^sxn^s)'/" 


y ^ ds(l - b^sin^e)^'^^ 


[K(b^) - E(b,)j 


(D-29) 


where the complete elliptic integrals are defined by Eq. (D-29) . They are 
well-studied and well tabulated functions. 

We have thus expressed the geometric mean beam length at the center 
of the side of a cylinder in terns of its radius R and length L. The trans- 
parent flux follows from Eqs. (D-9) , and an approximation for the flux over 
the whole range from black to transparent is given by Eq. (D-11) . 

This expression can new be used to find the at any point on the 
side of the cylinder. First ve recognize that if we wish to find the flux 
at the end of the side of the cylinder, rather than the nid-point of the side, 
that flux IS the same as half the flux from the mid-point of a cylinder twice 
as long. We merely replace L by 2L at^i 
D-29) . The result is 



j .D5 


divide by 2 in Eqs. (D-27, D-29 and 



,r /-TT/S 

I = ^ / d8 
0 IT y 


sin^S 


j^(L/2Rl^ + cos^ej^''^ 


2L [K{b) - ECb)] , 


= [(L/2R1^ + l] 


- 1/2 


(D-30) 


Now consider any point a distance x from the mid-point of the cylinder. 
The flux is made up of contributions from the end of the side of two cylinders, 
one of length L/2 + x and the other of length L/2 - x. 


i = - 
0 TT 


(§-)/ 


TT/2 


dS- 


2 

sin B 


[(L/2 + x)V4R^ + cos^b] 


1/2 


+ ^ 


,T ^■^/S 


dB- 


n 

sin' 6 


[(L /2 - x)^/4R^ + cos^B 


1/2 


(D-31) 


This can be used to express the flux at any point in terms of elliptic 
integrals, by writing it as 


V • I ! '— V * --- ['=<‘■.1 - 

b^ = [(L/2 ± x)^/4R^ + ij 


(L/2 - X) 


[K(b_) - E(b_)j| 


1 - 1/2 


(D-32) 


The resulting transparent flux is given by Eg. (D-9) , and the approximate 
flux at any by Eq. (D-11) . 

To find the power radiated out the wnole side of the cylinder, whose 
area is 2 ttRL, we use Eq. (D-12) , where dA = 2TiRdx. We thus need the integral 

_L/2 

^dx 

-L/2 


^0 = ^ 


Ij/ ^ 


of £- from Eq. (D--31) . If we let y = L/2 + x in the first integral and 
° i ( . • 

y = L/2 - X in md sedond integral, they both become the same, and we find > ■ 
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4 ^TT/ 2 2 /* ^ 

-±J dSsinBy 


y dy 

2 2 1/2 
(y/2Rl + cos 3 

. 


If we now let z = fy/2R)^, dz = y dy/2P , then 

,2 ^/2 


T , 8R 
0 “ TIL 


y* d 3 sin^S j 


L^/4R^ 


dz 


Q l^Z + COS^sj 


1/2 


The z integral is simple, and we find 


I = M 

0 ttL 


T TI/2 

- J dB sin* 3 


.1/2 


2 2 2 \ 

dB sin‘B ||l /‘IR + cos t j - cosBj 


The second term integrates to 1/3. The first can also be expressed as com- 

2 2 

plote elliptic integrals, by replacing cos 3 by 1 - sin ? and introducing 
the parameter b already defined in flq- (D-30) . The integral is 


f -1 2f 22l 
J 1*3 b Ein 2l 1 - b sin 3 


= G(b)/3 


J3(b) = (1 - b^) K(b) - (1 - 2b“) E(b) 


(D-33) 


The final form for the area-averaged geomotric mean beam length for the side 
of the cylinder is then 

1 - 1/2 


- _ IGR 

"^o ~ 3 ;il 


-1 
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G(b) , 
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ItJ * > 

^3 ■ ^ 
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(D-34) 


Again, the transparent radi'atcd power from the side is found from 
Eq. (D-14) as 


‘RcJB 


= 2-RLk 3 i 
U (J 0 


(D-35) 
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while the black radiated power is from Eq. (D-13) , 


Q„ „ = 2TrRLB 
Rofl ti) 


The approximate expression, Eq. (D-15) , becomes 

-k I 


Q„ = 21TRLD 
Rai (i) 


/, WO \ 

[1-. ) 


(D-36) 


(D-37) 


Eq. 


For the total surface of the cylinder, we car. easily find 
(D-19) as 



from 


= 2R/(1 + B/L) 

tot 


(D-38) 


since the total area is 2 ttRL + 27R'“. 
radiated powers are 


'‘RufT. 


tot 


The corresponding transparent and black 

47TR^Lk B 
0) (U 


A 


and 


"ROB 


tot 


= 27TRL(1 + B/L)B 

(J 


while the approximate form (D-15) is, with (D-38) 



= 27rRL(l + R/L)B 

L. 



(D-39) 


The use of Eq. (D-37) is still sorewhat inconvenient because of the 
elliptic integrals in the expression for G contained in 4^. Since the use 
of Eq. (D-37) is only approximate, some further simplification is 3 ustified. 

We note that the ratio of the total to side black power is gust the area ratio 
1 + R/ L, -while for the transparent power tnis ratio is 2R/Z^. A calcu^lation 


4 r 

_0 _ G(b) 

2R ” 37R ^3 

L D 


1 



• ( 


of 



shows that for R/L 1, this function is within 5.2% of (1 + B/L) always 
lying above it. So a satisfactory approximation is 



2R 

1 + B/L 


(D-40) 


Thus foi cylinders longer than their radius, the ratio of total to side 
transparent radiation is also approximately the area ratio. (For > R/L <_ 1 
the error is loss than 3.2%.) This implies that the sideways power is much 
larger than the base power for such cylinders. 


For R/L > 1 , we can expand the elliptic integrals for b near 1, 
and find the scries for G in the form 



G = 0,75(1 

- [' 

1 - |(i - v] 

In [l 6/(1 - K 2 )] 


+ 1 - 

2.25(1 ■ 

- K^) + 0.25(1 

- K^)*- 

4m 


K, - = 1 

(L/2R)^ 

-r 


This expression 

is accurate 

to 1%. 

(For 2R/L > 1 

It IS accurate 


lD-41) 


We conclude that the expressions to be used in the simple models of 
laser heated thrusters described in Appendices A and B are Eq. (D-39) for 
the power radiated from all surfaces of the cylinder, and Eq. (D-37) for the 
power radiated from the side surface, with given by Eq. (D-40) or (D-41). 
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APPENDIX E 


FLOW EQUATIONS FOR THRUST CHPiffiER DESIGN COMPUTER PROGRAM 


The laser-heated thrust chamber is represented by the axisymmetric 
flow flow of a gas which is in chemical equilibtium, is absorbing laser 
radiation, and has energy transport by convection, diffusion and radiation. 

A set of equations for this flow was presented in Section 5.1 of Ref. E-1. 

An alternate, and perhaps more convenient version of the energy equation will 
be presented here, and contrasted with the previous formulation. 

We consider the flow in an axisymmetric channel, whore axial grad- 
ients are small compared to radial gradients. This enables us to make the 
boundary layer approximation that the viscous transport terms need only in- 
clude radial gradients, not axial ones. Let us denote the axial and radial 
coordinates by x, r and the corresponding velocity components by u, v. The 
mass conservation equation is then 


9 pu ^ ^ 9rpv 
jx r 9r 


(E-la) 


The boundary layer approximation applied to the radial momentum equation 
shows that the pressure is constant over a cross-section, so that p = p(x). 
Then the axial nomenttun conservation equation is 


pu 


9u 

9x 


+ 




r 



(E-lb) 


The energy conservation includes the convection of enthalpy, the radi- 
ation, the flow work, the viscous dissipation and the conduction and diffusion 
of enerci’. Tf h denotes enthalpy, and ^ the radiation energy flux vector, then 
the conservation of energy requires ' 
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+ 




3h ^ ^ „ 

P” ^ ^ P'' 37 + 


dx 




r 3r 


^cf If ■" 21pv*i If^ 


(E-lc) 


The conduction term involves the frozen conductivity which includes the 
translational and internal contributions, but not the reaction contribution. 

It is the sura of the terms given in Eqs. (4-55a) and (4-55b) of Ref. E-1. 

The diffusion term is the summation involving the enthalpy per unit mass h^ 
of the i-t!i species, the diffusion coefficient for that species, and the 

gradient of the corresponding mass fraction C^. This form of the diffusion 
term is the simplified one usually used, where the species diffusion velocity 
has been expressed by Kick's law as 3C^/3r, to avoid the considerable com- 
plications of multicomponent diffusion. 

There are two ways to express the conduction and diffusion terms in 
terms of enthalpy. We start with the definition of the mixture enthalpy as 

h = Eh C 
1 1 


and differentiate to get 


dh = Eh dC + c ,dT, c , = Ec . C , c . = dh /dT 
1 1 pf pf pi 1 pi 1 


(E-2) 


Here we have recognized that for each species, depends only cn T, and its 
derivative is the species specu fic heat while c^^ is the frozen specific 

heat of the mixture, being just the mass-weighted sum of the specific heats of 
the individua I components . 
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One way to express the conduction-diffusion heat flux Lor.n to sub- 
stitute from Eq. (E-2) for dT, and combine the dC^ term with the diffusion 
term. The result is 


, 3T ^ u ^Cf 

A_- T — + JIpD.h. — = 

Cf 3r 1 i ^Tr c - 

pf 


c - L3r 

r\r ^ 


+ 1(1x3,-!) h, 


where the species Lewis number Loj^ is deiincd as 


LCi = PDi 

This approach is suitable for chemical non-equi librium flow, where the are 
computed from chemical reaction equations. It also shawr the simplification 
which occurs when Le^^ = 1 and the mass fraction gradients vanish. This form 
of the conduction-diffusion term was used in the flow equations written in 
Section 5.1 of Ref. E-1. 

A second form of the conduction-diffusion teji. can be used for chem- 
ical equilibrium flow, when is a function of p and T. Then we may write 
Eq. (E-2) IS 

dh = ^ Zh^ (dCj^/dT)p + o^.-j dT + Zhj^ (3Cj^/dp)^dp 

The first simmation on the right is the reaction specific heat and the sum 
of the terms in square brackets is the equilibrium specific heat 


c = c , + Zh (DC /3 t) 
pe pf lip 


Then the enthalpy derivative is 


dh = c vlT + Zh (DC /Dp) do 

pe 1 1 T * > 


We similarly use 


dC = (DC /Dt) dT + (DC /3p)_ dp 

^ X O X X 


(E-4) 

(E-5) 




f 


in the diffusion term, so the conduction-diffusion tern becomes 


''Cf 


3t 

3r 


ZpD^h^ 


3r 



Epo.h. 
1 i 




If ^ 


lE 

3r 


(E-6) 


The quantity in square brackets on the right is the equilibrium thermal 
conductivity 


(%), 

which adds the reaction conductivity term to the frozen conductivity. The 
reaction conductivity is the expression in Eq. (4-55c) of Ref. E-1, and the 
discussed there is the equilibrium conductivity X^^. If we now replace 
the temperature gradient in Eq. (E-6) by the enthalpy gradient from Eq. (E-5) , 
we find 


i ^ 

Cf 3r 


IpD.h. 
i 1 


3C. 

1 

3r 


X 


Ce 


c 

pe 


3h 

3r 


+ 




3r 


« 


The second term on the right involves an equilibrium Lewis number 


Le .V = pD c /X 
ei 1 pe Ce 


but It also has a radial pressure gradient term, which we have taken to be 

small enough to ignore. Thus, by using the equilibrium formulation of X^ 

and c , we can write the conduction-diffusion term simply as 
P 
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'cf 3r 


EpD^h^ 


3r 


Ce 3r 


Ce 


pe 


3r 


The energy conservation equation (E-lc) now becomes 

\2 , 


pu 


liL 

3x 


f~57 ♦ ’-S "“5S+>\37j*?37 (,'r;5?j 


The laser energy flux term is split off from the radiation flux vector S. 

We take a parallel beam with only an x component, which we call I. Then the 

laser provides a term 3l/3x on the left of the energy equation. However, the 

absorption equation for the laser energy is expressed in terms of a laser 

absorption coefficient k as 

L 


3l/3x = -k I 

Li 


This shows that the laser can be looked upon as providing a heating term k I 

L 

on the right side of the energy equation. (We permit I to depend on both r 
and X, through the laser propagates only in the x direction.) The remaining 
part of the radiant energy flux is due to gas radiation transport, and we 
will denote it by ^ . Then the energy equation is 


^ 3h . „ _K ^ /BuV 


+ k.I + 
L 


1 ^ 3h \ 

r 8r dr) 


(E-8) 


This is the equilibrium form of the energy equation. If the first 
approach to the convection-diffusion term had been used, Eq. (E-3) , the only 
difference would have been replacement of the last term by 


1 _a_ 

r 3r 


Cf 

'pf 


r ^ 

[3r 


+ Z (Le^^ -X) h 


3C n 

— ^ 
i 3r Jl 


(E-9) 


This form was given in Eq. (5-ld) of Ref. E-1. 
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l£ Le^ ■■ 1 for all species, the equations become the same , because 
it is easy to see that then *» ^e^^pe' equilibrium form 

Eq. (E-8) avoids the calculation of Sc^/Sr when the assumption Le^ « 1 is 
not warranted. This simplification is achieved at the expense of having to 
use instead of X^^, which means including the reaction conductivity. 

The present computer program uses the equilibrium formulation Eq. (E-8) , 
and the subsequent transformations described in this appendix will use that 
equation also. The alternative formulation using Eq. (E-9) is given in 
chapter V of Ref. E-1. 

These conservation equations are singulau: on the axis r ° 0, and must 
be expressed in limiting form there. On the axis, symmetry requires that 
V = 0 and 3u/3r = 0. Near the axis these two quantities are proportional to 
r. All the terms with the operator v3/3r vanish because v vanishes while 
3/3r is finite. The only terms which must be looked at carefully are those 
with a 1/r factor. A typical one is the viscous term in the momentum equa- 
tion (E-lb) : 


1 _ 3 _ 
r 3r 


rvi 


3 “ 

3r 



+ 


V ^ ^ ^ 

r 3r 3r 3r 


■ 2 

The last term is of order r and vanishes on the axis. The first term is 
finite, and the second term is the same as the first when a Taylor series 
«»xpansion of 3u/3r is made. Thus on the axis 


/ 1 3 

3u\ 

( 3^u \ 


\ r 3r 


V 3^2 j 



O 

a 


Ebcactly similar 

treatment 

can be accorded the second term in Eq. 

(E-la) and 

the last term in Eq. (E-8) 

. The resulting equations on the axis 

are 

3pu 

3x 

- 

0 

(E-lOa) 
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(E-lOb) 


3u 

3x 




dx 


3 2 
r 


^ 3h ^ „._R dp ^ T ^ « 3^ 

pu -s— + v-s =U-;^+k,I + 2 r- 

3x — dxL C.2 

pti 3r 


(E-10c>“ 


The numerical solution o£ the flow equations Is accomplished by using 
the stream function as an Independent variable. We define it as 


'I' 


I Our dr, 

Jo 


3ili 3tl/ 

-37= ^ 


■pvr 


so the mass conservation equation is satisfied automatically. The stream 
function has the physical significeince of being the mass flow through a 
cross-section of radius r, except for a factor of 2 tt. It replaces the vari- 
able r. The other independent variable isC= x, but is given ' new name 
to avoid confusion in the transformation process. The transformation from 
x,r to is given in Section 5.2 of Ref. E-1. The momentum and energy 
conservation equations (E-lb and E-8) become 


3u 

35 


. Jl§S. + 
pu dC 



2 

pur y 


3u \ 

3ii»; 


(E-lla) 


3h ^ V .S _ ^ ^ ^ 

35 pu p d5 pu 


2 


2 

pur p 



+ 


_3 

3i^ 



Ce ^ 
d\i> 


) 


(E-llb) 


The corresponding equations on the axis are found from Eqs. (E-lOb and 

E-lOc) . The transformation on the axis is derived in Section 5.2 of Ref. E-1 
2 2 

and shows that 3 /3r ■* pu3/3i|i. Then the momentum and energy conservation 
equations on the eucis are 


3u 

35 


1 dp 
pu d5 


+ 2p 


3u 

3\|) 


(E-12a) 
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9u 

ac 


+ 


(E-12b) 


7-_g ^ 

pu ” p d5 pu dt{i 


Tlie only other change to the flow equations given in Chapter V of 
Ref. E-1 is in calculating the heat flux to the wall eurising from the con- 
vective and diffusive terms. The radial heat flux is the coefficient of r 
in the last term of Eq. (E-lc) : 


q = - rpD.h - 

Cf 3r 11 


3C 


_1 

3r 


The minus sign is needed so that negative radial gradients produce positive 
radial flux. The transformation to the equilibrium form leads to the form 
exhibited by the last term of Eq. (E-8) : 


''Ce 3h 

~ c 3r 
pe 


(E-13) 


The further transformation to tp as the independent variable produces 


- pur 





(E-14) 


as in the last term of Eq. (E-llb) . It is shown in Section 5.6 of Ref. E-1 

how to evaluate this on the wall streamline where u = 0 but d/dip is infinite. 

2 2 
We convert 3/3r in Eq. (E-13) to 3/3r by multiplying by 2r. Then 3/3r 

is evaluated in finite difference form and related to a finite difference 

form in The result is 


2r X 
w Cew 


c 

pew 



eH_ 

4(}p -ip) 
w 


(E-15) 


where the quantities with subscript w are to be evaluated on the wall Stream- 
line and those without subscripts are to be evaluated at the first streamline 
away from the wall. 


as 



< 

We note that if the wall streamline is at an angle to the x-axis, i 

only the con^xjnent of the radial flux which is normal to the wall supplies I 

heat to the wall. | 

In sunarary, the computer program solves Eqs. (E-lla and E-llb) , , 

I 

except on the axis where Eqs. (E-12a and E-12b) are used. The conduction- 

diffusion heat flux to the wall is evaluated using Eq. (E-15) . The change 

from the formulation given in Chapter V of Ref. E-1 is in the energy equation , 

and q„ where X„,/c , is replaced by X /c , and the 3c./9r term is omitted. 

C Cf pf Ce pe 1 

The equilibrium thermal conductivity is defined in Eq. (E-7) , and its 
calculation is discussed in Section 4.4 of Ref. E-1. The equilibrium spec- 
ific heat Cp^ is defined in Eqs. (E-2 and E-4) . It is actually (3h/9T)^ 
for a mixture which stays in chemical equilibrium, as shown by Eq. (E-5) . 
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APPENDIX F 


RADIATION MODEL FOR MIXTURES 

A band model was constructed In the Interim Report, Ref, F-1, for the 
radiation from a specific mixture of hydrogen, water and cesium. The model 
has been modified to give reliable predictions of the radiation from mixtures 
subject to the following constraints; temperature ^ 6000 K; pressure between 
three and one hundred atmospheres; radius between three and thirty centimeters; 
initial cesium molar concentration from .005 to .01; and initial water molar 
concentrations from .05 to .10. 

The radiating species which were considered were H^O, OH, Cs (bound- 
free transitions!, (vibrational transitions), H (bound-free transition), H , 
cind the free-free contributions from E-H, E-Cs, E-H^ and E-Cs^. Analysis of 
detailed calculations of both the absorption coefficient and the radiation 
emission coefficient indicated that the three bands defined in Ref. F-1 were 
adequate to describe the radiation for the range of parameters considered. 

The parameters used to describe the radiation in the three bands are described 
below, 

4 -1 

Band 3:2x10 cm <a)<®* 

The major contributions in this band are the free bound radiation 
involving H and the first excited state of Cs. The expression used in 
Ref. F-1 is still valid, that is 

k* = l.OE-17 n* + 2.5E-17 n„. (F-1) 

3 Cs HH 

where n* is the number density of the first two excited states of cesium 
Cs 

and n,_ is the numbe?- density of H . 

KM 

’ -14-1 

Band 2; 8000 cm < o) < 2 x 10 cm (F-2) 

The major contributions in this band are the free bound radiation 
involving H and the intermediate excited states of Cs. The expressions 
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developed in Ref. F-I remains valid, that is 

k* = 1.75E-17 n- + 4.2E^17 

where n** is the number density of the intermediate excited states of cesium. 
Cs 

Band 1: 0 < (0 < 8000 cm"^ 


There are several important radiators in the low temperature band. 

In addition to bound-free radiation associated with highly excited states 
of Cs, there are contributions from Bremsstrahlung when an electron collides 
with Cs, H, Cs^ and H^, and also contributions from vibrational bands 
of OH and H^O. The contribution of the cesium bound-free radiation and 
the Bremsstrahlung radiation are well represented by the following expres- 
sion for the absorption coefficient 

k* = 1.3E-17 n*** + r n n^O^.to = 4500 cm"^) (F-3) 

wL Cs { JL E 


*** 

In the above expression n^^ is the number density of the highly excited states 
of cesium, n^ is the electron concentration, n^ is the concentration of species i 
and IS the inverse Bremsstrahlung cross-section for collision between 
electrons and species i evaluated at (o = 4500 cm 


The absorption coefficient for the vibrational transitions in H^O, 

OH and H^ varies rapidly as the wave number varies; thus these transitions 
are not amenable to accurate representation by a constant absorption coefficient. 
Fortunately, for the range of radius, pressure, temperature and composition 
considered here, these transitions can be well represented as a transparent 
contribution to the radiation. They are included in the radiation modeling 
as a contribution to V*S , namely 




1.07E-36 T°‘^ 

604 3/T 
e - 1 


2.53E-11 

5136/T , "oH 
e -1 


(F-4) 
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where n n , and n are the number densities of HO, H_ and OH 
H O OH 2 2 

respectively. 

The band model calculations are compared to detailed calculations 
of the radiation parameters in Figs, F-1 - F-6, In Fig. F~1 the spectral 
absorption coefficient for 1 atm of Cs at 6000 K is shown for both models. 
The band model values are average values of throughout limited spectral 
ranges; however, the average values depart from the detailed values both 
at very low (j and at U) greater than 30000 cm”^. In Fig. F-2 the normalized 
spectral emission coefficient is shown in both models. It is seen that the 
major contribution to the emission coefficient is well approximated, and the 
errors in k^ at both small and large 0) do not lead to errors in the spectrally 
integrated emission coefficient of more than 5 percent. 

The absorption coefficient for one atmosphere of H^ at 4000 K is 
shown in Fig. F-3 for both the detailed model and the band model. The band 
values are representative of average values except for (jJ > 30,000 cm ^ and 
to < 3,000 cm~^. The major contribution to the spectral emission coefficient 
is again well approximated; the error in the band model approximation for 
the integrated emission coefficient is 10 percent. 

Finally, Fig. F-5 shows the detailed calculation of the absorption 
coefficient for vibrational transitions in one atmosphere water vapor at 2000- K. 
The rapid variation of the absorption coefficient cannot be well characterized 
by a single constant value from O) = 0 to 8000 cm but the integral of the 
spectral emission coefficient, which is shown in Fig. F-6, is given by a term 
in Eq. F-4. This approximation is in error only when strong reabsorption 
occurs, which requires that the product of the partial pressure of H^O and 
the radius be approximately 30. Only for the largest sizes and highest pres- 
sures considered does the transparent approximation fail, and even then the 
failure is limited to low temperature regimes where radiation losses are small. 

Based on similar comparisons of the band model with detailed spectral 
calculations, it has been concluded that the band model approximates the im- 
portant radiative losses within ten percent for the range of pressure, temper- 
atuie, radius, and compositions of interest. 
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Figure F-2 Normalized spectral emission coefficient and 
band model approximation for cesium. 
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APPENDIX G 


LASER SUPPORTED COMBUSTION WAVE MODEL 

One method of coupling laser energy into a gas is by means of a laser 
supported combustion (LSC) wave. This is a mechanism which permits the ab- 
sorption of laser energy in a flowing gas, to produce a high temperature flow. 
It is analagous to a deflagration wave in chemical combustion. The LSC wave 
consists of a constant pressure region in which the gas temperature rises 
sharply to a high value, and then decays slowly. The high temperature region 
is maintained by the cibsorption of laser energy, usually by inverse Brems- 
strahlung absorption. The rapid temperature rise from the cold on-coming 
gas to the hot region is accomplished by the transport of energy from the 
hot region forward into the cold gas by both thermal conduction and radiation. 
The decay of the temperature at the back of the wave is caused by radiation 
losses dominating the absorption of laser energy, because the beam has already 
been almost completely absorbed. 

A previous analytical study has been made of such waves in hydrogen, 
and the results have been reported in Ref. G-1. Presented there is a pri- 
marily one-dimensional model of the structure of such a wave, although the 
effect of beam convergence and some estimate of transverse energy losses 
were included. 

In the present study, the node! in Ref. G-1 was partially extended 
to include an equilibrium mixture of cesium and hydrogen. These extensions 
were incorporated into the computer program developed in the previous work, 
so that the structure of I.SC waves in cesium-hydrogen mixtures can bo calcu- 
lated. This appendix gives a brief description of the model as it presently 
exists. 

For an LSC wave, as for a deflagration wave, the momentum equation 
is satisfied by taking the pressure p to be constant. Conservation of mass 
in the one-dimensional flow requires 
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Pu «• m, 
A 


« constant 


(G-1) 


where p and u cire the density and flow velocity. Conservation of energy is 
conveniently expressed by the use of the total enthalpy as 


m. 


_d_ 

dx 


(h 


+ u^2) 




(G-2) 


Here h is the gas enthalpy, I the laser intensity, the gas radiation 
energy flux, T the temperature, R the beam radius, and the equilibrium 
thermal conductivity of the gas, as defined in Appendix E. This equation 
expresses the balance of convection of total enthalpy, laser absorption, 
gas radiant energy flux, combined conduction-diffusion energy flux, and an 
estimate for tr 2 uisverse conduction losses. 

The relation for the beam intensity is obtained by conserving the 
power in the beam. If the beam area is then the power is lA^, and con- 
servation requires the decrease of power in the propagation direction x to 

equal the power absorbed with absorption coefficient k . Thus, we have 

x* 

• -dIAg/dx = kj^IAg (G-3) 


whose solution is 






k dx 


(G-4) 


where the subscript i refers to the location x = 0. 


The energy equation (G-2) deals with the flow near the beam centerline 
where dl/dx = -k I. Use of Eq. (G-4) shows that 


dl 

dx 


-k 


L 


-T 

e 



(G-5) 
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Thus, the laser absorption term depends not only on the absorption coefficient, 
but also on the beam shape (x) . 

In Ref. G-1, Section-4.2,-an-expre3slon-f or- the- beam-shape is obtained 
in terms of three physical parameters. One is the distance x^, from the be- 
ginning of the wave, x ° 0, to the beam focal plane. The other two are the 
dlstcince L from the focusing lens to the focal plane, and the size u) of the 

la 

beam on the focusing lens. Thev appear only in the ratio L/u , which is the 

Xj 

f- number of the lens. The area ratio is given by 


O F 


(G-6a) 


where 


■ ^ (i) 


(G-6b) 


and li) is the laser wavelength. Equations (G-5) and (G-6) serve to define 

Xi 

the laser absorption term in the energy equation (G-2) in terms of the 
absorption coefficient k^, -the wave to focal plane distance x^ and the lens 
f- number L/u . 

Xi 

The radiation term V*£ must be specified in terms of the gas prop- 
erties. The study for hydrogen radiation, made in Refs. G-1 and G-2, showed 
that there were two parts to the radiation transport. One was a -volumetric 
loss term consisting of transparent radiation, P^, and one a radiation trans- 
port term which can be modeled using the conduction approximation. Then 
the radiation term is expressed as 


VoS*^ = P_ - 


T dx 


T 

(vi)-^/ V 
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A parametric expression for was given in Appendix D of Ref. G 
for hydrogen as 


17.2 X 10 


.9 V 




A 

elA 




(G-8) 


3 

expressed in W/m . Here 6 is the ionization energy for hydrogen, expressed 

^ -3 

in degrees K, n, is the number density of hydrogen atoms in m and Q , 

A elA 

is the electronic partition function of these atoms 


?1A 


E 2 exp j -0 (l-j“^) /T 
j=l 


The constants and depend on gas pressure amd laser power level, and 
have been worked out for 10 kW aind 5MW lasers at 1,3,10 and 30 atm: 


P(kM) 

p(atm) 

S 

_f2 

10 

1,3 

2 

2.4 

10 

10,30 

1.3 

2.4 

5000 

1,3 

2 

0.4 

5000 

10,30 

2.2 

0.4 


For the radiation conductivity X , an expression was also given in 

R 

Appendix D of Ref. G-1. It is expressed for hydrogen in terms'of (in MKS 
units) 


X 


vuv 


0.775x10 


30 


0^ ,, ,-150,865/T 
"a (11605/T)® 


7 

^ (150,865/T) 


(G-9a) 




(G-9b) 
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Then it was shown that 


X <the smaller of 1 or 2Rk,„„,l 

R VUV \ VUVj 

is to be used in the derivative term of Eq. (G-7) and 
R VUV 


(G-lOa) 


(G-lOb) 


is to be used in the integral term. 

The constants C in Eq. (G-2) and C in Eq. (G-7) are coefficients 
LC X*R 

which provide the strength of the transverse conduction losses. This form 
of transverse loss estimate was first proposed by Raizcr in Ref. G-3, for 
ordinary conduction. He proposed the value 3 for the constant, which we have 
adapted for C . For C we have used 3 for a 5 MW laser and 0 for a 10 kW 
laser, where the size of the beam is quite small for the intensities of 
interest. 

With the laser and gas properties so specified, the energy equation 

(G-2) can be integrated from x=0, at which we choose a temperature low 

enough so no radiation losses cam occur up to that temperature. A boundary 

condition there is obtained by equating the energy which has been added to the 

incoming gas to that which has been transported from the hot gas. This energy 

balance requires, with ( ) the incoming gas state, 

o 

m, [h - h + (u ^ - u “)/2l = (X^ + X„) (dT/dx) (G-11) 

A ‘ 1 o 1 o Ce R 1 1 


The equation of state for the gas is also used, relating p, p, T and h. 
Finally, the mass conservation equation (G-1) permits u to be eliminated in 
favor of p. 

The integration of Eq. (G-2) requires a value for m^. This is a 
free parameter whose value must be chosen so that the integration produces 
a physically reasonable temperature profile. If m^ is too small, T will go 
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to zero after rising to a peak. If is too large, T will climb continually 
without reaching a peak. For the correct solution, T rises to a pe^lk, and 
then decays to a constant value. The procedure for finding the correct 
value of tn , by iterative integration of Eq. (G-2) , is described in Ref. G-1. 

A 

The formulation as described above is also valid for a gas mixture, 
provided the gas properties and radiation and conduction loss coefficients 
are properly evaluated. For the cesium-hydrogen mixture of present interest, 
some of these properties have been evaluated, but not all. The equilibrium 
thermodynamic properties have been described in Chapter IV of Ref. G-1 and 
are incorporated in the present program. The equilibrium thermal conductivity 
was also considered in Chapter IV of Ref. G-1. There it was found that small 
amounts of cesium added to hydrogen (Cg/H=0.01) have a negligible effect on 
the thermal conductivity of the mixture, so that the values of appropriate 

for pure hydrogen may be used. They are given in Fig. 4-1 of Ref. G-1, 

J ^T 

X dT. 

300 ^ 

The only part of the program which has not yet had the effect of 
cesium incorporated is that arising from radiation transport, both and 
X . The present version of the program still contains the formulation appro- 
priate to pure hydrogen. It~was found in Chapter III of Ref. G-1 that 
cesium is a good radiator, so it will be necessary to modify the radiation 
terms to include the presence of cesium seed. This has not yet been done. 
Therefore, the use of cesium in the present version (September 1980) of the 
LSC wave program must be viewed with great caution. The radiation transport 
cind loss terms will probably be under-estimated by an unknown amount. These 
radiation terms need to be re-formulated to maikc them valid for inclusion of 
cesium. 
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APPENDIX H 

MODEL FOR LASER SUPPORTED COMBUSTION 
WAVE HEATING WITH BUFFER GAS COOLING 

One of the possible heating inodes for a laser heated thruster is use 
of the laser supported coisbustion (LSC) wave mechanism. In this mode nearly 
all the laser energy absorption occurs in a few centimeters in an LSC wave, 
and a hot gas issues from the wave, with very little energy remaining in the 
laser beam. This hot gas must then be accelerated to produce thrust. 

In a pure hydrogen LSC wave, the ten 5 >erature of the hot gas is near 
20,000 K, while for a cesium-seeded gas it may be somewhat lower. In either 
case, the gas is too hot to be in contact with the walls. One possible cool- 
ing method is to surround the hot gas with a co-flowing unheated annulus of 
buffer gas, which will act as an insulator between the laser-heated gas and the 
wall. This method of cooling can be modeled by combining the LSC wave model 
described in Appendix G with the thrust chctmber design model described in 
Appendix E. 

A diagram of the model is shown in Fig. H-1. An LSC wave of radius 

Rj^ stands at the entrance to the thrust chamber. Entering the wave is a laser 

beam of intensity I^ and a mass flow per unit area m^ = P^^u^ caused by cin 

incoming stream of gas at temperature T^, density p^, pressure p^ and speed 

u . Also entering the thrust chamber is buffer gas flow, with the same prop- 
o 

erties and speed, in the annulus from R^ to R^. However, this gas is not 
heated by the laser, since the laser beam is only of radius R^. 

At some point in the structure of the LSC wave we arrive at station 
A-A where we stop the LSC wave calculation and start the two-dimensional thrust 
chamber calculations. At present, this s'lation is taken to be at the peak tem- 
perature of the LSC wave, though other choices could be made. It is desirable 
to have nearly all the laser energy absorbed by station A-A, or else further 
increases in ten 5 >erature could occur in the thrust chamber. 
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The conditions at the station A-A in the LSC wave (velocity, tem- 
perature, pressure) are used as initial conditons for the start of the two- 
dimensional thrust chamber design calculation in the radius region up to 
In the annulus between and R^, the initial conditions are-determined, from 
the buffer gas flow. However, the incoming buffer gas conditions themselves 
are not used. Instead, the energy lost by transverse conduction from the LSC 
wave to station A-A is added to the incoming buffer gas flow. The resulting 
gas conditions are used as initial conditions in the thrust chamber design cal- 
culation in the annulus from R. to R • 

n m 

The energy lost, per unit thickness of the LSC wave, in the transverse 
direction, is modeled in Appendix G as 

T T 


^ ^Lc r 

^"7 J 




/ 




(H-1) 


per unit transverse surface area. 


Thus, the total energy lost is 
A-A 


Q = TTR 


/ 


5a 

dx 


dx 


(H-2) 


where the integral is taken- up to station A-A. 

This energy is assumed to be absorbed by the buffer gas, which has a 

2 2 

ma^s flow rate of m^ per unit area, through an annulus of area ~ Rj^) • 

Then the additional energy per unit mass added to the buffer gas is 




= Ah 


(H-3) 


An energy balance for the buffer gas from inlet to station A-A is 


then 


h 

o 


+ u^/2 + ah = hj^^ + u^j/2 


(H-4) 
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where the subscript bf refers to the buffer gas at station 
a mass balance requires 

*A = Po% = ^bf'^f • 

The pressure is assumed not to change through the LSC wave 

the caloric equation of state "^bf^ thermal 

Pl^(P t T. ,) are assumed known for the buffer gas mixture, 
bt o or 

H-5) are sufficient to determine the state at station A-A when the incoming 

state ( ) is known, 
o 

This procedure corqiletely specifies the state of the gas at station 
A-A, the initial station of the thrust chamber design calculation. It con- 
sists of a hot core and a cool buffer annulus, both made of the same gas 
mixture. Then the flow equations specified in Appendix E can be applied to 
calculate the thrust chamber shape and heating. 

It should be noted that the present (September 1980) version of the 
thrust chamber design computer code does not contain a radiation model for 
the high temperature (> 6000 K) gas emerging from the LSC wave. It should 
therefore be used with caution until such a model has been included. 


A-A. Similarly, 

(H-5) 

so p, , = p . Since 
bt , o 

equation of state 
Eqs. (H-3, H-4 and 
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APPENDIX I 


EXPERIMENTAL EMEASUREMENTS OF ABSORPTION OF 
10.6 pm RADIATION BY WATER VAPOR 

- Ebcperimental Methodology ^md limitations 

An existing shock tube facility was used to measure the absorption 
coefficient of water vapor for 10.6 pm radiation at high temperatures. The 
shock tube is a well understood experimental tool and is ideal for creating 
the pressures and temperatures required for the proposed measurements. 

Small shock tubes are relatively inexpensive to operate and take the gas 
mixture almost instantly to a known high temperature and pressure. Test 
times up to a millisecond in duration are possible, which are more than 
sufficient to simulate the heating process in a laser-heated thruster. 

As shown in Fig. I-l, a shock tube consists of a high pressure sec- 
tion and a low pressure section separated by a diaphragm. When the diaphragm 
bursts, the high pressure gas flowing into the low pressure section acts as 
a piston, compressing and accelerating the low pressure gas ahead of it. Pro- 
cessing of the low pressure gas is accomplished by a shock wave which travels 
ahead of the piston. When the shpek wave passes through the test section at 
the end of the shock tube, it reflects off the end wall and reshocks the pre- 
viously shocked gas, compressing it further. It is usual to specify the con- 
ditions in the high pressure region (the driver gas) by the subscript 4, 
conditions in the undisturbed low pressure section (the driven gas) by sub- 
script 1, the gas conditions behind the incident shock wave by subscript 2, 
and behind the reflected shock wave by stabscript 5. This is shown in the 
position versus time plot (Fig. I-lb) . Measurements are made in the test 
section in regions labeled 2 and 5. 

The shock wa/e heats the processed (or driven) gas, resulting in an 
enthalpy per unit mass Ah, which is easily and accurately calculcible when 
the shock velocity is known. The enthalpy of the processed gas is approx- 
imately doubled when tne shock wave is reflected off the end wall. High 
temperatures require a high enthalpy per unit mass, thus implying high shock 
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THE SIMPLE SHOCK TUBE 




(a) Conctnttonal ihock tube (b) An (x—t) diagram showing progress of the 
shock Mare, the rarefaction fan end the contact surface separating drirer and experi- 
mental gases (c) The pressure distribution along the tube at time /, (d) The temperature 
distribution at tune /, 


Fig. I-l (from Ref. I-l) 




volocitloa. Tho velocity of the gas piston increases as the pressure ratio 

across the diaphragm increases but Ccuinot exceed a value of a few times the 

2 

initial sound speed in the driver gas a, whore y is the specific 

heat ratio and a is the sound speed in the driver gas ) , even with an infinite 
pressure ratio. Therefore, driver gases with high sound speeds are desired. 
Most shock tube drivers use hydrogen, and in the PSI shock tube, room temper- 
ature hydrogen at up to 200 atmospheres pressure is used. 

An idealized representation of the temperature, pressure and laser 
energy transmitted which will bo observed at the tost section is shown in 
Fig. 1-2. Until t^, when the incident shock* passes, the low values of 
temperature and pressure of the quiescent driven gas are observed, together 
with the incident laser Intensity I^. When the incident shock passes, at t^, 
the temperature and pressure rapidly rise to the higher values and p^. 

The transmitted laser intensity, I^, is reduced bccauoe of the absorption 
in the shocked gas. When the reflected shock passes, at t^, a second jump 
in temperature and pressure occurs, to T. and p. . The transmitted laser in- 
tensity again drops because of further absorption in the hotter, higher 
pressure gas behind the reflected shock. 

When performing an experiment in water vapor, there are limitations 
on the choice of initial pressure in the driven gas (region 1 in Fig. I-lb) . 

The initial pressure in this low pressure section of the shock tube is de- 
termined by the allowable range for the partial pressure of water vapor. 

The low pressure limit is determined from the requirement that any condensed 
water on the shock tube walls must not affect the measured water vapor mole 
percent of the gas in the tube. At low partial pressures of water vapor 
there will be at most a monolayer of water due to chemisorption, and wo have 
estimated (using theories of Langmuir, Brunauer, Emmett and Toller in Ref. 1-2) 
that van der Waals forces also yield a monolayer of water molecules on the 
shock tube walls. To eliminate any effect of this monolayer on the gas in 
the tube, we have calculated that the partial pressure of water vapor in the 
driven gas must bo greater than 10 ^ atm. 
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t. - arrival time of incident shock 
‘ at test station 


t - arrival time of reflected 
^ shock at test station 


t 


1 profiles of temperature, pressure, and 
intensity for test station near end wall. 



similarly, if the pressure in the low pressure section is raised 
to the point vdiere the vapor pressure of water exceeds the saturation vapor 
pressure at the coldest point in the system, water will condense out at that 

point, again making it difficult to control the amount of water in the tube. 

-2 

Therefore, the maximum vapor pressure of water is limited to 2 x 10 atm, 
the saturation vapor pressure at 20°C. 

Within these high and low pressure limitations, we have determined 
the range of temperatures and pressures attainable in the shock tube by per- 
forming a series of calculations for the water vapor mixtures using chemical 
equilibrium computer codes to determine the pressure and temperature behind 
the^incident and reflected shocks (regions 2 and 5 in Fig. I-lb) . The re- 
sults of these calculations are presented in Figs. 1-3 and 1-4. Calculations 
were made of the pressure and temperature behind the incident and reflected 
shocks for numerous water vapor mixtures. The mole percent of the water 
vapor was varied from 2% to 16%. An example result for 10% H^O in is 
presented in Fig. 1-3. In this plot the pressure P 2 denotes pressure be- 
hind the incident shock and p^ denotes the pressure behind the reflected 
shock. The pressure and temperature behind these shocks are increased by 
increasing the driver pressure (and thus, the shock strength) . The upper 
limit to the pressures is reached with the driver pressure at 200 atm. Lines 
are drawn corresponding to the two limiting values of the initial total pres- 
sure in the low pressure section (i.e., p = .01 atm corresponds to the 

-3 ^ 

10 atm H^O partial pressure minimum) . By connecting the maximum pressure 

points for the two limits on p^^, we can define operating regimes for the shock 

tube. These regimes are shaded in Fig. 1-3 for both the incident and reflected 

shocksi 

It can be seen that the mixture of H^O and permits operation up 
to 1600 K behind the incident shock and up to 2900 K behind the reflected shock. 

A similar plot for a mixture of 10% H^O in argon is shown in Fig. 1-4. 
This shows a much increased temperature range, up to more than 6000 K. It is 
easier to raise the temperature of the heavy gas argon than the light hydrogen 
because the speed of sound is lower in argon. 
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Fig, 1-3 pressure-tenperaturc region which can be reached 
in a 10% H^O/90%H2 gas mixture. 
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Fig. 1-4 Pressure temperature region which can be 
reached in a 10% 11^0/90% Ar gas nuNture. 
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The preliminary experiments reported here were performed in water- 
argon mixtures, to take advantage of the higher temperatures which could be 
reached. For laser propulsion applications, hydrogen is the most desirable 
propellant so water-hydrogen mixtures cire of the most interest. This brings 
upthe question of the relation between the absorption coefficients measured 
in water-argon mixtures and those which would occur in water-hydrogen mix- 
tures at the same temperature and with the same partial pressure of water 
vapor. 

At low temperature and atmospheric pressure the absorption coefficient 
IS strongly dependent upon the location of tlie laser line being absorbed 
relative to the center of the molecular absorption lines of the water vapor 
molecule. The width of these absorption lines are dependent on the collision 
partner, i.e., the diluent gas. This would seem to indicate that the measured 
cibsorption coefficient could be strongly dependent upon whether argon or 
hydrogen is the collision partner. However, the number of populated states, 
and the density and the widths of the molecular absorption lines increase sub- 
stantially at elevated temperatures and pressures. Therefore, the lines 
overlap and the collision partner dependence of the absorption coefficient 
will become less and less important as the pressure and temperature is in- 
creased. 

Under the present program, experiments were made using the P(20) laser 
line, with argon as the collision partner. The experiment is being continued 
at PSI under different sponsorship, and the dependence of water vapor absorption 
coefficient on both collision partner and laser line will be explored in that 
program. 

The shock tube operating regime must be combined with the constraints 
resulting from the minimum measurable amount of the absorption which takes 
place. This constraint places an additional limit on the range of allowable 
water vapor pressures. There must be sufficient absorption by the H^O vapxir 
so that it can be accurately measured. The minimum absorption we can measure 
is 1% of the probe laser intensity. For the shock tube diameter of 3.8 cm, 
we then need an absorption coefficient k such that 3.8k = 0.01. 
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From the literature, we have estimated the absorption coefficient 
of water vapor (Ref. I- 3 ) and the results are shown in Fiq. 1-5. From 
this figure we can find the partial pressure of water vapor at any tempera- 
ture which will give a value k = 0.01/3.8 = 0.00263 cm Figure 1-6 shows 
the result of such a calculation. Two curves are shown; one assumes that 
the water molecules remain water molecules even at high temperatures; (that 
is, the chemistry is "frozen") and the other assumes that the chemistry is 
in equilibrium. The latter curve should apply everywhere except very near 
the shock front. With this curve the required vapor pressure starts low 
at room temperature because the water vapor continuum is helping absorb, 
and rises until lines become important; then the required pressure falls. 

It continues to fall until the temperature reaches 2500 K, at which point 
the water molecule begins to dissociate and the required initial water 
vapor pressure must rise again. For frozen chemistry the pressure also 
rises slightly at high temperature. By operating the shock tube with 
partial pressiires of water above the equilibrium curve in Fig. 1-5, more 
than 1% cibsorption will occur, which should provide sufficient measurement 
accuracy. 

Absorption by water molecules condensed on the windows, through 

which the absorption measurement will be made, places another limit on 

the permissible operating partial pressures of water. Liquid water has a 

large absorption coefficient for 10.6 )im radiation, with a value of k 
3 -1 

greater than 10 cm ; therefore, water must not bo allowed to condense 
onto the windows. Although the initial gas sample will have a sufficiently 
low partial pressure of water, the shock wave raises this pressure by large 
factors. This is, of course, no problem in the middle of the shock tube 
whore the gas is hot; however, the gas near the wall, which is at this same 
high pressure, is at or near the temperature of the wall. Since the satura- 
tion vapor pressure of water is a steep function of temperature, warning the 
test section wall can prevent condensation. The vapor pressure of water as 
a function of temperature is as follows: 




Vapor pressure p^ = 0.02 atm 0 293 K 
0.12 atm 0 323 K 
1.0 atm 0 373 K 
2.5 atm 0 400 K 
4.7 atm 0 423 K 

From these values of water vapor pressure, we observed that if the 
wall temperature could be maintained around 400 K, then condensation would 
not occur near the windows under experimental conditions. This approach was 
used early in the program; however, the shock tube overheated, resulting in 
the melting of several piezoelectric pressure transducers and O-ring vacuum 
-seals. Programmatic constraints did not allow for a major heating system 
redesign so the heating was reduced to 310 K. At this temperature, the par- 
tial pressure of water vapor behind the reflected shock is much greater than 
the saturation pressure and the effect of condensation (in the cold boundary 
layer near the windows) had to bo investigated. The results of the worst 
case conditions are presented below. 

In the shock tube the separation distance between the shock and the 
contact surface (region 2) approaches a limiting value because of the wall 
boundary layer (Ref. 1-4). In addition, the cold boundary layer, if it is 
thick enough, will absorb laser radiation in a manner different from the 
uniform hot gas in the center of the tube. It is necessary to estimate the 
contribution of this Ixiundary layer to the overall absorption measurement. 

Studies of the wall boundary layer in a shock tube have been re- 
ported by many investigators (Refs. 1-5, 1-6 and 1-7). Here we will use the 
Mirels and King series expansion solution (Ref. 1-7) to estimate the wall 
boundary layer contribution to the absorption of laser radiation. Consider 
the case of a free sLieum temperature of approximately 600 K in region 2, 
and a pressure of one atmosphere. This can bo considered the most serious 
case for boundary layer contributions because the absorption coefficient of 
water vapor at the laser wavelengt); has a minimum at about 600 K (Fig. 1-5) . 
Thus, the absorption coefficient in the boundary layer is higher than in the 
center of the tube. The mixture is taken to be 20% H^O and 60% H^. The 
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Fig. 1-6 Conditions which produce 1% cibsorption in water vapor. 
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Fig. 1-7 Shock tube boundary layer temperature 
profile at X = 4.2 cm. 
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c.<perimental test time will be ~ 20 Psec with a shock speed of approximately 
2.1 X 10^ cm/sec (from a shock tube calculation). We will calculate the wall 
boundary layer profiles at a distance of 4.2 cm from the shock. 

From the zeroth order solutions, we find that the boundary layer 

-2 

thickness at 4.2 cm behind the shock is approximately 4.5x10 cm (2.4% 
of the shock tube diameter) . In order to estimate the wall boundary layer 
contribution to the overall absorption of laser radiation, the temperature 
profile at X = 4.2 an, where X is distance behind the shock, is also needed 
and shown in Fig. 1-7. We can evaluate the integral n = /k dy where 
k(T) and y are the absorption coefficient of water vapor and the normal 
distance from the wall, respectively. For the free stream portion (in the 

-4 -5 

tube center) , n = 7. x 10 . The wall boundary layers give t) = 1.2 x 10 

and this is approximately a 1.7% effect. The above calculation is for a 
pressure of 1 atm. As the pressure is raised to 50 atm while the tempera- 
ture is maintained at 600 K, the wall boundary layer contribution to the 
laser radiation absorption becomes approximately 12 percent. Therefore, 
we conclude that boundary layer effects will not seriously influence the 
curacy of the absorption measurement. 

Experimental Apparatus 

figure 1-8 is a photograph of the 1.5" I.D. shock tube during the 
set-up period for the NASA water vapor absorption coefficient measurements. 

In the foreground is the driver, connected to the 22' long driven section. 

At the far end there is a 10" long test section where measurements are made 
and finally a 9 cubic foot dump tank is used to bleed the tube to atmospheric 
pressure after a shot. The line tunable OT CO^ laser and the optical de- 
tection system are located on the table perpendicular to the tube at the 
test section. Behind the teible is the power supply for the laser. The laser 
absorption measurement system will detect a 1% change in absorption. (The 
heating tapes to raise the wall temperature had not been installed when the 
photograph was taken.) Also located on the tiible but not readily visible 
is a 0.25 m Jarrell-Ash Spectrometer used for laser wavelength calibrations 
and a Scientech power meter to measure the laser output. 
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^^9: The PSI 1,5 Inch Shock Tube Facility 


Figure Z-9 1& a schematic representation of tbs' shock tube gas 
handling system. Gas mixtures used in a test are prepared in the mixing 
tank. The tank Is heated to prevent condensation of H^O, and can bo evac- 
uated to a hard vacuum with the diffusion dump. The water is added from a 
flask through a valve at one end of the tank. The water vapor pressure is 
monitored' by a Validyne DP-7 variable reluctance pressure transducer. 
Diluent gas is' then added to the water vapor by injection through a spray 
bar located on the central £ucis on the tank, and the total pressure of 
the gas is measured. 

The pressure in the driven section of the shock tube is monitored 
by a thermocouple gauge before filling. Gas from the mixing tank is valved 
into the tube and the pressure is monitored by the Validyne transducer. 

The driver is evacuated to a rough vacuum before pressurization. 
Driver gas is supplied by a six cylinder manifold, and the fill pressure 
is monitored by a Bourdon-type pressure gauge equipped with, a check valve 
to hold the diaphragm burst pressure. When the diaphragm bursts, the shock 
wave is created and propagates through the driver section toward the tost 
section. The pressure behind the shock wave was originally measured by 
seven piezoelectric transducers at 2.5 ft intervals upstream of the tost 
section, and the shock velocity was determined by measuring the time inter- 
val between successive pressure transducers. After the tube overheated 
and damaged the transducers, the shock speed was determined by measuring 
the shock transit time between one transducer upstream of the test section 
and the main transducer at the tost section. 

The test section and absorption measurement system are located at 
the end of the driver section and arc shot/n schematically in Fig. I-IO. 

The output of a PSI-built line-tuned single mode OT CO^ laser is directed 
through the test section to measure the absorption in the shocked gas. The 
windows in the test sections were 3/4" in diameter by 0.2" thick zinc sul- 
fide. Originally NaCl windows were used; however, the thermal shock caused 
surface crazing after only one shot. 
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riq. 1-9 Schematic of shock tube gas handling system 





A HgCdTe detector was used to monitor the absorption of the hot gas 
during a shot. Before the shot, the beam was blocked to obtain the 0% 
transmission level; no absorption is observed in the cold gas. A Kistler 
piezoelectric pressure transducer measures the pressure history in the test 
section during the absorption measurement. An end wall is located 1.5" 
downstream from the windows to create the reflected shock; this allowed 
sufficient time for measurements behind the incident shock. Idealized 
traces of the tem'perature , pressure and absorption in the test section are 
shown in Fig. I-IO. 

After the experiment, the gas in the tube automatically bleeds into 
the dump tank through a check valve, eliminating the need for a downstream 
diaphragm. The data was collected by a CAMAC data acquisition system 
and processed by the PSI PRIMB 400 computer after each shot. 

Experimental Results and Discussion 

Shock tube runs were conducted in H^O/Ar mixtures to measure the 
absorption coefficient of water vapor. In the early checkout runs, large 
fluctuations in the transmitted laser beam intensity were observed, with a 
10 to 20 Ms period, after the gas was shock-heated. 

Two possible sources of these fluctuations were identified. One 
was a 7% increase in the cross-sectional area of the shock tidse in the test 
section, where flat windows were mounted in the walls of a round tube. The 
passage of the incident and reflected shock waves through this area generates 
weak secondary waves which cculd cause small density changes in the hot gas. 
This in turn could cause changes in the index of refraction of the gas, move 
the transmitted laser beam around on the detector, and cause fluctuations 
in the signal measured. 

A second possible cause could be a parasitic coupling of the distur- 
bances in the hot gas to the gas in the laser cavity, causing fluctuations 
in the laser output. 
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Fig. I-IO Schematic diagram of the IR absorption measurement system 
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In the present program It was not possible to design and build a new 
test section to eliminate the area change, so Both problems were attacked by 
a reconfiguration of the optical system. The collimated laser beam was dif- 
fused by inserting a roughened NaCl flat in the beam before it entered the 
test section. A lens was used to image the center of the test section onto 
the detector. The optical diffuser decouples the laser from the test gas to 
reduce the possibility of parasitics being generated. Diffusion of the laser 
beam illuminates the test section more uniformly and reduces the effect of 
local density fluctuations on the beam. 

This optical configuration reduced the fluctuations to a few percent 
of the transmitted beam intensity. However, it also reduced the transmitted 
intensity sufficiently that broad-band emission from the hot gas in the test 
section, after the reflected shock, was comparable to the original laser in- 
tensity. To reduce this emission, the optical bandwidth of the detection 
system was reduced by inserting a 10.6 Pm (1.7 pm FWHM) interference filter 
in front of the detector. The hot gas emission was substantially reduced, 
but not eliminated by the filter. To further reduce emission effects, the 
data runs were made by performing two identical shots at each experimental 
condition, one with the laser off to measure gas emission, and one with the 
laser on to measure both gas emission and transmitted bean intensity. Sub- 
traction of the detector signal for the two runs yielded the actual -absorption. 

The temperature and pressure of the test gas was calculated using 
the measurement of the incident shock velocity. This measurement was made 
in the uniform portion of the shock tube upstream of the tost section. 

However, the gas in which the absorption was observed was processed after 
the shock moved into the test section, where the area increased by 7%, which 
caused the shock to slow down. A small correction was applied to the 
measured shock velocity to account for this slowdown. 

The correction was obtained from the theory of shock dynamics in 
channels of varying area (Ref. 1-8) , which expresses the change in shock Hach 
number from an original value to a new value M, when the area changes from 
an original value A^ to a new value A. The expression for the new Mach number 
is 
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The temperature and pressure changes are then found by perturbing the equa- 
tions for temf>erature and pressure ^ump across a shock wave. 

A typical computer output trace of the pressure measured in the test 
section is shown in Fig. I-ll, with the jumps at the incident and reflected 
shocks clearly visible. Also shown is the detected signal with the probe 
laser off. The signal shows a slight increase when the reflected shock 
arrives, due to hot gas emission. The second shot of this pair of shots, 
with the laser on, is shown in Fig. 1-12. The pressure trace is nearly 
identical. The 100% laser transmission line is shown, together with the 
measured trace. A small amount of reduced transmission (absorption) is 
seen after the incident shock, and a much larger amount after theTreflected 
shock. The figure shows that the level of the transmitted signal still 
varies, so a value was determined by averaging over 50 Ps of digitized datfi 
(50 points in Fig. 1-12) . 

Five shots have been made by the method described, and for each shot 
a measurement was made behind the incident shock and the reflected shock. 

The results are shown in Fig. 1-13. The points labeled 2 are incident shock 
measuicments, and those labeled 5 are reflected shock measurements. These 
results, and their relation to the values used in the modeling, are discussed 
in the main body of this report in Chapter IV. Table I-l presents the re- 
sults in tabular form. 
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Fig. I-ll Experimental pressure and detector traces with probe laser off. 
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TABLE I-l 

EXPERIMENTAL PARAMETERS 


Run 

Gas Con 5 >osition 



A 

16.7% HjO in Ar 

B 

1.00% HjO in Aj* 

C 

4.96% HjO in Ar 

D,E 

10.3% H^O in A].* 


Incident Shock Conditions 


Run 

T^ (K) 

(atm) 

Ph2o 

» Trans 

k (cm 'ama '^) 







A 

1148 

1.69 

0.0664 

>99 

<0.040 

B 

1106 

11.3 

0.0566 

>99 

<0.047 

C 

1435 

7.02 

0.0663 

>99 

<0.032 

D 

1721 

4.50 

0.0733 

97.2 

0.10 

E 

1681 

4.38 

0.0729 

98.0 • 

0.073 


Reflected Shock Conditions 


Run 

Tg (K) 

Pg (atm) 

Ph2o 

% Trans 

k (cm ^ama ^) 

A 

2101 

8.04 

0.171 

94.2 

0.092 

B 

2146 

49.0 

0.123 

93.7 

0.14 

C 

2770 

33.4 

0.141 

78.5 

0.45 

D 

3064 

22.8 

0.162 

77.9 

0.41 

E 

2999 

22.1 

0.167 

77.8 

0.38 
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The error bars shown are those due to an estimated 2% error in 
reading the transmitted laser intensity. They do not account for signal 
fluctuations, or possible errors in the temperature and pressure correc- 
tions duo to the area change in the tost section. 

As mentioned above, this experiment is continuing at PSI under 
other sponsorship. In that work a new test section will be used, and other 
experimental problems will be eliminated. The present results should be 
looked upon as preliminary. They do, however, provide encouraging evidence 
that the absorption coefficients used in the modeling are within a factor 
of 2 of the actual values. 
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